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We derive the equations of time-independent stochastic quantization, without refer- 
ence to an unphysical 5th time, from the principle of gauge equivalence. It asserts that 
probability distributions P that give the same expectation values for gauge-invariant ob- 
servables (W) = J dA W P are physically indistiguishable. This method escapes the 
Gribov critique. We derive an exact system of equations that closely resembles the Dyson- 
Schwinger equations of Faddeev-Popov theory. The system is truncated, and solved non- 
perturbatively, by means of a power law Ansatz, for the critical exponents that characterize 
the asymptotic form at k « of the gluon propagator in Landau gauge. For the tranverse 
and longitudinal parts, we find respectively D T ~ (k 2 )~ 1 ~ otT ~ (A; 2 ) 043 , suppressed and 
in fact vanishing, though weakly, and D L ~ a (k 2 )~ 1 ~ OLL ~ a (/c 2 ) -1 ' 521 , enhanced, with 
cut = —loth- Although the longitudinal part vanishes with the gauge parameter a in the 
Landau-gauge limit, a — > 0, there are vertices of order a -1 so, counter-intuitively, the lon- 
gitudinal part of the gluon propagator does contribute in internal lines in Landau gauge, 
replacing the ghost that occurs in Faddeev-Popov theory. We compare our results with 
the corresponding results in Faddeev-Popov theory. 
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1. Introduction 

1.1. Some recent developments in non-perturbative QCD 

The problem of the strong interaction presents an exciting challenge. One would like to 
understand how and why QCD describes a world of color-neutral hadrons with a mass gap, 
even though it appears perturbatively to be a theory of unconfined and massless gluons 
and quarks. Clearly an understanding of non-Abelian gauge theory at the non-perturbative 
level is required. Happily, there has recently developed a convergence of results by dif- 
ferent methods: (i) non-perturbative solutions of the truncated Dyson-Schwinger (DS) 
equations in Faddeev-Popov theory, (ii) numerical evaluation of gauge-fixed, lattice QCD 
propagators, and (iii) exact analytic results. The agreement between these very different 
methods almost 5 decades after the appearance of the original article of Yang and Mills [1] , 
would indicate that by (ii) we are beginning to get reliable values of the gluon propagator 
in the unbroken phase, and by (i) an understanding of the mechanism that determines 
it. This motivates the present investigation in which we derive the DS equations of time- 
independent stochastic quantization and solve them by truncation and a power-law Ansatz 
for the gluon propagator in the asymptotic, low-momentum regime. In accordance with 
earlier results by methods (i), (ii), and (iii), we find that, compared to the free propa- 
gator l//c 2 , the would-be physical, transverse component of the gluon propagator is short 
range, while the unphysical, longitudinal component is long range. 

As concerns (i), solutions of the DS equations, the decisive step was taken in [2], 
where a solution of the truncated DS equations in Faddeev-Popov quantization in Landau 
gauge was obtained for which the transverse gluon propagator is short range, while the 
ghost propagator is long range. These properties were confirmed in subsequent DS calcu- 
lations, using a variety of approximations for the vertex [3], [4], [5], and [6]. More recent 
calculations extend the asymptotic infrared and ultraviolet solutions to finite momentum 
k, without angular approximation [7], [8]. All these calculations 1 give a transverse gluon 
propagator in Landau gauge D T (k) that is highly suppressed in the infrared compared 
to the free massless propagator l//c 2 , and that in fact vanishes lim/j^o D (k) = 0, at 
k = 0, in some cases weakly, like a small positive power of k. Indeed, according to the 
present calculation it vanishes like (/c 2 ) 043 . A review of DS equations in QCD may be 

1 Stingl [9] had earlier obtained a solution of the DS equation with the property that D T (k) 
vanishes at k = 0, without however including the ghost loop, whereas the ghost loop gives the 
dominant contribution in the infrared region in the recent solutions. 
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found in [10]. In the present work we shall discover a close conection between the ghost 
propagator in Faddeev-Popov theory and the longitudinal part of the gluon propagator in 
time-independent stochastic quantization. 

Concerning (ii), numerical studies, it is striking that an accumulation of numerical 
evaluations of the gluon propagator in Landau gauge also show qualitive suppression of the 
gluon propagator at low momentum, both in 3-dimensions on relatively large lattices, [11], 
[12], [13], and in 4 dimensions, [14], [15], [16], [17]. Suppression of the gluon propagator and 
enhancement of the ghost propagator at low momentum has been reported by [18], [19], and 
[20]. Similar numerical results were obtained in Coulomb gauge, where an extrapolation 
to infinite lattice volume of the 3-dimensionally transverse, would-be physical, equal-time 
gluon propagator Dij(k) was consistent with its vanishing at k = 0, [21]. In QCD in 
the Coulomb gauge, the instantaneous Coulomb propagator, -D44, is closely related to the 
ghost or Faddeev-Popov propagator, and is a strong candidate for a confining potential. 
Significantly, D44 was found to be long range [21]. 

A recent numerical calculation in the Landau gauge, [16], reports a finite value of 
D T (k) at k = 0. This is strongly suppressed compared to l//c 2 , and suffices to exclude a 
free massless gluon. It might be thought that the finite value of D T (0) reported in [16] 
contradicts the zero value, D T (0) = 0, found here. However it is difficult to distinguish 
numerically between a finite value at k = and one that vanishes weakly, like (/c 2 ) 7 , with 
a small value for the infrared anomalous dimension such as 7 = 0.043 found here. For 
this function is almost constant down to very low k, and then veers toward zero with an 
infinite slope. Moreover a numerical determination of the continuum propagator at k = 
requires an extrapolation to infinite lattice volume. To establish a discrepancy it would be 
necessary to take 7 as a fitting parameter, and determine the numerical uncertainty in this 
quantity after extrapolation to infinite lattice volume, and this has not been done. Present 
numerical and analytic results are not inconsistent, within the considerable uncertainty of 
the numerical extrapolation to infinite lattice volume, and both agree that there is strong 
suppression compared to 1/k 2 . 

Infrared suppression of the gluon propagator D(k) and enhancement of the ghost 
propagator G(k) in Landau gauge was first found by Gribov, using avowedly rough ap- 
proximations [22]. He obtained the formulas, D{k) = k 2 /[(k 2 ) 2 + M 4 ], and, in the in- 
frared, G(k) ~ l/(/c 2 ) 2 , by restricting the region of functional integration to the interior 
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of the Gribov horizon in order to avoid Gribov copies. 2 In Coulomb gauge he also ob- 
tained a long-range Coulomb potential. Concerning (iii), exact analytic results, it was 
subsequently found, [25] and [26], that restriction to the interior of the Gribov horizon, 
enforced by a horizon condition, yields at k = 0, both the vanishing of the gluon propa- 
gator limfc^o D(k) = in Landau and Coulomb gauge, and the enhancement of the ghost 
propagator lim^o k 2 G(k) = oo. 3 

It was at first surprising that the solution of the DS equations obtained in [2], [3], [4], 
and [5] agreed with these exact results that are a consequence of cutting off the functional 
integral at the Gribov horizon, for this condition was not imposed in solving the DS 
equations. However it was subsequently pointed out [6] that the DS equations in Faddeev- 
Popov theory depend only on the integrand, and the fact that the integral of a derivative 
vanishes provided only that the integrand vanishes on the boundary. The key point is 
that the integrand does vanish on the Gribov horizon for the Faddeev-Popov determinent, 
det[— D(A) ■ d], vanishes there (as explained in footnote 2). Thus Gribov's prescription to 
cut off the functional integral at the (first) Gribov horizon, is not a constraint that changes 
the DS equations, but rather it resolves an ambiguity in the solution of these equations [6]. 
The cut-off at the first Gribov horizon assures that both the gluon and ghost Euclidean 
propagators are positive, which is a property of the solutions obtained for the truncated 
DS equations. Moreover the solution of the DS equations in Faddeev-Popov theory with a 
cut-off at the Gribov horizon is the only one for which a comparison with numerical gauge 
fixing to the lattice Landau gauge is (approximately) justified. For as explained in footnote 

2 We remind the reader that numerical gauge-fixing to the Landau gauge is achieved by mini- 
mizing, with respect to local gauge transformations g(x), a lattice analog of -Fa(<7) = / c? 4 x| 9 A| 2 . 
At any minimum, this functional is (a) stationary, and (b) the matrix of second derivatives 
is non-negative. These conditions correspond to (a) the Landau gauge condition d ■ A = 0, 
and (b) the positivity of the Faddeev-Popov operator —D(A) ■ d which, moreover is symmetric 
—D(A) -d = —d-D(A), for d- A = 0. Condition (b) defines the Gribov region, so numerical studies 
of the Landau gauge automatically select configurations within the Gribov region. Positivity of 
— d ■ D(A) means that all its eigenvalues A n are positive, and the boundary of the Gribov region, 
known as the (first) Gribov horizon, is where the first (non-trivial) eigenvalue vanishes. Thus the 
Faddeev-Popov determinent, det[— D(A) ■ d] = Yl n ^ n > which is the product of the eigenvalues, 
is positive inside the Gribov horizon and vanishes on it. These considerations do not apply to 
numerical gauge fixing to the Laplacian gauge [23] , [24] . 

3 It is noteworthy that the confinement criterion of Kugo and Ojima [27] and [28] also entails 
lmifc^o k 2 G(k) = oo. 
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2, numerical gauge fixing to the Landau gauge automatically produces a configuration that 
lies inside the Gribov horizon. Thus a consistent picture emerges of the gluon and ghost 
propagators in QCD using the different methods (i), (ii) and (iii). 

1.2. Diffulties of Faddeev- Popov method at non-perturbative level 

The DS calculations [2] [5] rely on Faddeev-Popov theory which however is sub- 
ject to the well-known critiques of Gribov [22] and Singer [29]. At the perturbative level, 
Faddeev-Popov theory is unexceptionable, and elegant BRST proofs are available of per- 
turbative renormalizability and perturbative unitarity [30] . In lattice gauge theory however 
the BRST method fails because the total number of Gribov copies is even, but they con- 
tribute with opposite signs, leading to an exact cancellation [31], [32]. In continuum gauge 
theory, the Faddeev-Popov-BRST method may nevertheless be formally correct at the 
non-perturbative level without a cut-off at the Gribov horizon, if ones sums over all signed 
Gribov copies [33] , [34] . However even if this is true, it would imply large cancellations be- 
tween copies, that may amplify the error of an approximate non-perturbative calculation, 
and even the Euclidean gluon propagator D(k) is not necessarily positive. Alternatively, 
one may choose the solution of the DS equations in Faddeev-Popov theory that corre- 
sponds to a cut-off at the first Gribov horizon, which indeed is our interpretation of the 
solutions of [2] - [8]. Hopefully, this is an excellent approximation. But it remains an ad 
hoc prescription that is not correct in principle because of the existence of Gribov copies 
inside the Gribov horizon [35] and [36]. 

Wilson's lattice gauge theory provides a quantization that is both theoretically sound 
and well suited to numerical simulation. It also provides a simple analytic model of con- 
finement in QCD by giving an area law for Wilson loops in the strong-coupling limit. A 
striking feature of lattice gauge theory is that both the numerical simulations and the 
strong-coupling expansion are manifestly gauge invariant. This manifest gauge invariance 
provides a paradigm for continuing efforts to understand confinement in QCD. Neverthe- 
less it may be worthwile to pursue other approaches. The vexing problem of bound states 
in quantum field theory is particularly urgent in QCD where confinement causes all phys- 
ical particles to be bound states of the fundamental quark and gluon constituents. In this 
regard it is noteworthy that even the simplest of all bound-state problems, the Hydrogen 
atom, is not easily solved in a gauge-invariant formulation. 
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1.3. Review of stochastic quantization 

In order to avoid the difficulties just mentioned of the Faddeev-Popov method, we 
turn to stochastic quantization of gauge fields for, as we shall see, this method provides a 
correct continuum quantization at the non-perturbative level. Stochastic quantization has 
been developed by a number of authors [37], [38], who have expressed the solution as a 
functional integral [39], and demonstrated the renormalizability of this approach [40], [41]. 
A systematic development is presented in [42], [43], [44], [45], [46], [47], reviewed in [48], 
that includes the 4-and 5-dimensional Dyson-Schwinger equation for the quantum effec- 
tive action, an extension of the method to gravity, and gauge-invariant regularization by 
smoothing in the 5th time. Renormalizability has also been established by an elaboration 
of BRST techniques [49], [50]. Stochastic quantization may be and has been exactly sim- 
ulated numerically including on rather large lattices, of volume (48) 4 , [51], [52], [53], [54], 
[55]. This suggests the possibility of a promising interplay of DS and numerical methods. 

In its original formulation [37], stochastic quantization relies on the observation 
that the formal Euclidean proabability distribution Pq(A) = iVexp[— Sym(A)], with 4- 
dimensional Euclidean Yang-Mills action Sym(A), is the equilibrium distribution of the 
stochastic process defined by the equation, 



^= f ( fix-*—(-¥-+*?™-P) (11) 
8t J 6A°(x)\6A°(x) 6A°(x) )' K > 

Indeed it is obvious that Pq{A) is a time-independent solution of this equation. Here t 
is an artificial 5th time that is a continuum analog of the number of sweeps in a Monte 
Carlo simulation of the Euclidean theory defined by the action Sym(A). As explained in 
sec. 3, this equation has the form of the diffusion equation with "drift force" — jj^f^j , and 
is known as the Fokker-Planck equation. The same stochastic process may equivalently be 
represented by the Langevin equation 

dA^ = JS Y M a , , 

dt 5A- [i ' Z) 

where A a ^ = A®(x,t) depends on the artificial 5th time, and corresponds in a Monte-Carlo 
simulation to the configuration on the lattice with points x M , with [i = 1, ...4 at sweep t. 
Here rf^ = rj®(x, t) is Gaussian white noise defined by (r]®(x, t)) = and (r]^(x, t)rj^(y, t')) = 
25(x — y)S flu 5 ab 5(t — t')). If A^exp[— Sym(A)] were a normalizable probability distribution 
- which it is not — every normalized solution to (1.1) would relax to it as equilibrium 
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distribution. However the process denned by (1.1) or (1.2) does not provide a restoring 
force in gauge orbit directions, so probability escapes to infinity along the gauge orbits, and 
as a result P(A, t) does not relax to a well-defined limitiing distribution lim^oo P(A, t) ^ 
N exp[— Sym(A)] . Nevertheless, according to [37], expectation values (0{A)) t of gauge- 
invariant quantities O(A) calculated at fixed but finite time t according to either of the 
above equations do relax to the desired Euclidean expectation value, (O) = linit^ 00 ((9)t. 

Unfortunately the renormalization program cannot be carried out in this scheme as 
stated, because that requires that gauge-non-invariant correlators also be well defined. A 
remedy is provided by the observation [38] that the Langevin equation may be modified 
by the addition of an infinitesimal gauge transformation, D a ^v c = {d^5 ac + f abc A b ^)v c , 



dAl 5S 
dt 



— + D?v° + rt. (1.3) 



Clearly this cannot alter the expecation-value of gauge-invariant quantities. Symmetry 
and power-counting arguments determine v a = a~ 1 d\A < ^ = a~ x d • A a , where a is a free 
parameter. For a > 0, the new term, that is tangent to the gauge orbit, provides a restoring 
force along gauge orbit directions, so gauge-non-invariant correlators also exist. 4 The new 
scheme is renormalizable. Only a harmless gauge-transformation has been introduced, 
so the Gribov problem of globally correct gauge-fixing is by-passed, and a continuum 
quantization of gauge fields that is correct at the non-perturbative level has been achieved. 
The modified Langevin equation is equivalent to the modified Fokker-Planck equation 

f -/ d4l 4w(^- A '^ )p )' (L4) 

where the "drift force" now includes the infinitesimal gauge transformation [38] , 

K >^-jA^) +a ~ lD ^ d - AC{x) - ^ 



4 To establish that the new force is globally restoring, we note that the hilbert norm of A is 
decreasing under the flow defined by the new force alone, A M = a~ l D^d- A. We have d\\A\\ 2 /dt = 
2{A fl ,A fl ) = 2a-\A f ,,D fl d ■ A) = 2 a - 1 (A^,d l ,d ■ A) = -2a,- l \\d ■ A\\ 2 < 0. This also shows 
that the region of equilibrium under this force is the set of transverse configurations, d ■ A = 0. 
Similarly, from d\\d- A\\ 2 /dt = 2(d ■ A, d ■ A) = 2a" 1 (d ■ A, d ■ Dd ■ A) it follows that this equilibrium 
is stable inside the Gribov horizon, where —d ■ D is a positive operator, and unstable outside it. 
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The additional "force" is not conservative, and cannot be written, like the first term, as 
the gradient of some 4-dimensional gauge-fixing action, a~ 1 D® c d ■ A c (x) ^ — 5 ^. s ^ ■ With 
this term, the normalized solutions P(A, t) to (1.4) do relax to an equilibrium distribution 
lim^oo P(A, t) — P(A), and Euclidean expectation values are given by the 4-dimensional 
functional integral, (O) = J dA 0(A)P(A). Although we cannot write P(A) explicitly 
because the force is not conservative, we do know that it is the normalized solution of the 
time-independent Fokker-Planck equation 



This equation defines what we call "time-independent stochastic quantization", and Hpp 
is called the "Fokker-Planck hamiltonian" . The solution P(A) of this equation provides a 
satisfactory non-perturbative quantization of gauge fields. 

[To avoid possible confusion of terminology, we note that stochastic quantization, 
whether in the time-dependent or time-independent formulation, — where "time" is the 
artificial 5th time — increases the number of dimensions by one as compared to the cor- 
responding standard Faddeev-Popov formulation of gauge field theory. Thus the solution 
of the time-dependent Fokker-Planck equation (1.4) can be usefully represented [39] as a 
functional integral with a local 5- dimensional action I = f dtd 4 x£$, whereas in Faddeev- 
Popov theory, expectation values may be calculated by a functional integral with a local 
4-dimensional action S = f d 4 x[(l/4)F^ u + ...]. Likewise the Fokker-Planck "hamilto- 
nian" Hpp determines, by the time-independent Fokker-Planck equation HppP = 0, a 
Euclidean probability distribution P(A) whose argument is a field A(x) that is a func- 
tion in 4-dimensional space-time with points x^, \i = 1, ...4. By comparison the quantum 
mechanical hamiltonian Hqm in ordinary quantum field theory determines, by the time- 
independent Schrodinger equation Hq M ^ = E^> , a wave-functional ^{A), whose argument 
is a field A(x) that is a function in ordinary 3-space x = (^1,^2,^3). Thus Hpp is not 
a quantum mechanical hamiltonian at all, but rather, it claims the name "hamiltonian" 
as the generator of time translations in the time-dependent Fokker-Planck equation (1.4), 
where the "time" is the artificial 5th time. Unlike the quantum-mechanical hamiltonian for- 
mulation, time-independent stochastic quantization is 4-dimensionally Lorentz (Euclidean) 
covariant.] 




(1.6) 
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Despite the development of stochastic quantization in [37] [50] it has apparently 
not so far been used for non-perturbative calculations in QCD, apart from [6]. 5 This may 
possibly be due to the complication caused by the extra "time" variable. Although the 
time-dependent formulation allows an elegant representation, with a local 5-dimensional 
action, it has the complication in practice that the gluon propagator depends on two 
invariants D(k 2 ,uj) instead of only one D(k 2 ). This prevents a simple power Ansatz for 
the infrared behavior l/(k 2 ) 1+a that allows one to determine the infrared behavior of 
the 4-dimensional theory self-consistently. For this reason we turn to time-independent 
stochastic quantization, where the correlators have the same number of invariants as in 
Faddeev-Popov theory. 

1.4- Outline of the present article 

We shall not use the 5-dimensional formulation here, but only the 4-dimensional, time- 
independent Fokker-Planck equation (1.6). The solution P(A) to this equation cannot be 
represented as a functional integral over a local 4-dimensional action. Nor shall we attempt 
to construct an explicit solution to (1.6). Our strategy instead will be to convert it into a 
system of tractable DS equations for the correlators. 

As a first step, we convert (1.6) into the DS equation, (6.3) below, for the quantum 
effective action V. The DS equation for T appears relatively complicated, with a second- 
order structure inherited from the second-order operator in (1.6). The main methodological 
innovation of the present approach is that the second-order equation for F is replaced, in 
sees. 6 and 7, by the much simpler DS equation (6.6) for a quantity, Q^(x), that we call 
"the quantum effective drift force". Indeed the new equation Q£(x) = K£(x) + (loop 
integrals), where K^(x) = —jjt^ + a~ 1 (D fl d ■ A) c , has the same structure as the first- 
order DS equation for V in Faddeev-Popov theory, j^F^ = sa^(x) + (1°°P integrals). In 
both of these equations, the leading term may be interpreted as a drift force and, most 
helpfully for the renormalization program, it is local in A(x). 

In the present work we give an improved treatment, as compared to [6], of the longi- 
tudinal degrees of freedom in the Landau-gauge limit a — > 0. In that work we integrated 
out the longitudinal degrees of freedom in the Landau-gauge limit a — > 0. This gave a 
time-independent Fokker-Planck equation for the transverse degrees of freedom only, with 

5 The equations of stochastic quantization have however been applied to dissipative problems 
in QCD, where t is the physical time, and x physical 3-space, [56]. 
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an effective drift force that was however non-polynomial and non-local. 6 By contrast, in 
the present work, the difficulty of a non-polynomial drift force is avoided by retaining the 
longitudinal degrees of freedom. Of course the longitudinal part of the propagator vanishes 
with the gauge-parameter a in the Landau-gauge limit lim — > 0. However the drift force 
(1.5) gives a vertex that diverges like 1/a and so, counter-intuitively, the longitudinal part 
of the propagator in the Landau gauge limit gives a finite contribution in internal loops, 
somewhat like the ghost in Faddeev-Popov theory. 

We shall be satisfied here to calculate only the infrared asymptotic form of the propa- 
gator, because that is where the challenging, non-perturbative confining phenomena man- 
ifest themselves. At high momentum, QCD is perturbative, and it has been verified to 
one-loop order by various methods [57]and [46], including the background field method [58], 
that stochastic quantization yields the standard /3-function. We leave for another occasion 
a numerical calculation which would be necessary to connect the high- and low-momentum 
limits. 

Since we use only the time-independent formulation here, we present, in sees. 2 and 3, 
a new derivation of eq. (1.6) that does not refer to the unphysical 5th time. At the end 
of sec. 3 the Minkowskian form of time-independent stochastic quantization is presented. 
[Some readers may prefer to go directly to sec. 4, which begins with (1.6).] The new 
derivation is more powerful, and yields new results, in particular, the Ward identity of 
Appendix C, and the proof in Appendix A that the kernel of the Fokker-Planck hamiltonian 
for quarks depends on gauge parameters only. We shall derive it from the obvious principle 
of gauge equivalence which asserts that probability distributions P(A) that give the same 
expectation values for gauge-invariant observables (W) = J dA W(A) P(A) are physically 
indistiguishable. We show that time-independent stochastic quantization provides a class of 
positive, normalized proabability distributions P(A, a), parametrized by a gauge parameter 
a that are gauge equivalent P(A 7 ai) ~ P(A 7 CI2), and that includes includes the Yang-Mills 
distribution iVexp(— 5Vm) as a limiting case. This method of quantization of gauge fields, 
in which the unphysical degrees of freedom are retained but controlled, is closely related 
to the physics of our solution of the DS equations. Indeed we find that the physical 
degrees of freedom are short range, whereas the unphysical degrees of freedom are not only 
present but of long range. In Appendix A, we extend the method to include quarks, and 

6 This was in turn decomposed into a conservative force that reproduced the Faddeev-Popov 
determinent, plus a second term that was neglected in the solution found in [6]. 



9 



in Appendix B, to lattice gauge theory. In Appendix C, we derive a Ward identity that 
controls the divergences of the theory. 

In sec. 8 we derive the explicit form of the DS equation for the gluon propagator. In 
sees. 9 - 11 we adopt a simple truncation scheme, and by means of a power-law Ansatz we 
solve for the infrared critical exponents that characterize the gluon propagator in Landau 
gauge asymptotically, at low momentum. The transverse part of the gluon propagator is 
short range, and the longitudinal part long range. In the concluding section we compare 
our results with calculations in Faddeev-Popov theory, and we interpret their qualitative 
features in a confinement scenario. We also suggest some challenging open problems, and 
possibilities for comparison with numerical simulation in lattice gauge theory. 

2. Gauge equivalence 

We first consider Euclidean gauge theory and later the Minkowskian case. Non-abelian 
gauge theories are described by the Yang-Mills action Sym(A) = (1/4) f d 4 x(F® ) 2 , where 
F£ v = d^A® ~d v A a ^ J rgf ahc A h [1 A c v . The Euclidean quantum field theory is formally defined 
by the probability distribution Pym(^4) = iVexp[— Sym(A)], and by the expectation- values 
(W) = J dA W(A) P Y m(A), normalized so (1) = 1. The challenge of quantizing a non- 
abelian gauge theory is that Pym(A) is not really normalizable because of the infinite 
volume of the local gauge group. 

The challenge would be hopeless, but for the fact that we are interested only in 
observables that are invariant under local gauge transformations, W( 9 A) = W(A), 
for all g(x), where 9 A fI = g~ 1 A fI g + g~ l d^g. This suggests the notion of gauge- 
equivalent probability distributions. Two probability distributions are gauge equivalent, 
Pi (A) ~ P2(A), if and only if (W)i = (W)2, for all gauge-invariant observables W, where 
(W)i = J dA W(A) Pi(A). Gauge equivalence of probability distributions is dual to 
gauge invariance of observables. Distributions that are gauge-equivalent are physically 
indistinguishable. Our solution to the quantization problem will be to replace the formal 
probability distribution iVexp[— Sym(A)], by a normalizable distribution that is gauge- 
equivalent to it. More precisely we shall exhibit a class of gauge-equivalent normalized 
distributions that includes iV~exp[— Sym(A)] as a limiting case. 
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3. A machine that makes gauge-equivalent probability distributions 

The construction of gauge-equivalent probability distributions relies on an equation 
that has the same form as the time-independent Fokker-Planck equation that is used to 
describe diffusion in the presence of a drift force. In this section, for simplicity, we deal 
with continuum gauge fields, or gluodynamics, only. The extension to quarks is given in 
Appendix A, and to lattice gauge theory in Appendix B. 

In order to simplify the appearance of various equations, we shall, as convenient, use 
the index notation A x , instead of A^(x), where the subscript x represents the triplet x, fi, a. 
We use discrete notation and the summation convention on the new index so, for example, 
f£ replaces J d 4 x . 

Let P(A) be a positive, P(A) > 0, normalized, J dA P(A) = 1, probability distribu- 
tion or concentration. In simple diffusion theory there is associated with this distribution 
a current, 

dP 

J x = -h— + K x P, (3.1) 

that is composed of a diffusive term, — fi-§j- : proportional to the gradient of the concen- 
tration, with diffusion constant h, and a drift term, K X P. Here K x is the drift force, as in 
Ohm's law with unit conductivity. We have introduced Ti for future convenience for a loop 
expansion which is an expansion in powers of Ti. Conservation of probability is expressed 
by the equation of continuity ^ = —§^-- The analogy of interest to us here is associated 
with the time-independent situation only. 7 In this case the current is divergenceless 

S = °' (a2) 

which reads 

H ™ p = ik ( - n (k + K * )p = - (3 ' 3) 

This is the time-independent diffusion equation with drift force K x . We call the linear 
operator defined here the Fokker-Planck "hamiltonian" , although Hpp is not hermitian, 



7 Stochastic quantization [37], including a drift force tangent to the gauge orbit [38], has 
traditionally been based on the time-dependent Fokker-Planck equation ^ = —HfpP, and relied 
on relaxation of the stochastic process to an equilibrium distribution that satisfies HfpP = 0. 
Here t is an additional, unphysical time variable that corresponds to computer time in a Monte 
Carlo simulation. By contrast, in the present article, the quantization of the non-Abelian gauge 
field follows from the geometrical principle of gauge equivalence, from which we derive the time- 
independent equation HfpP = directly, without reference to the additional time variable. 
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as it would be in quantum mechanics, and it certainly is not the quantum mechanical 
hamiltonian of the gauge field. 

We must be sure to choose a drift force that is restoring, so this equation determines a 
positive normalized distribution P(A). If the drift force were conservative K x = — 9 gJ M , 
then the normalized solution would be P(A) = iVexp(— Sym)- Gauge invariance of the 
Yang-Mills action, = 0, means however that the conservative drift force — g^f™) 

provides no restoring force in gauge-orbit directions. This is remedied by introducing an 
additional force K^ t jl (x) = D® c v c that is an infinitesimal gauge transformation, so the 
drift force is made of a conservative piece and a piece that is tangent to the gauge orbit, 

(3 - 4) 

^YM nac c 

5A«{x) + " ■ 

Geometrically, the drift force is a vector field or flow, and it is intuitively clear that a flow 
that is tangent to the gauge orbit has no effect on gauge-invariant observables. We will 
not fail to choose v a (x; A) so that D^v is a restoring force, to insure that (3.3) possesses a 
positive, normalized solution. Apart from this restoring property, v a (x; A) may in principle 
be an arbitrary functional of A. The time-independent Fokker-Planck equation reads 
explicitly 

H FP P= [d*x—?—-\- n-^+(-^L + D-v c )p] =0. (3.5) 

This equation is a machine that produces normalized probability distributions P V (A) that 
are gauge equivalent to iVexp(— Sym)- 

We now prove the basic result. Positive, normalized solutions of the diffusion equation 
(3.5) for different v are gauge equivalent P v ~ P v > , and include iVexp[— Sym (A)] as a 
limiting case. Our solution to the problem of quantizing a gauge field is to use any one 
of the P V (A) to calculate expectation- values of gauge-invariant observables. We consider 
observables that are invariant under infinitesimal local gauge transformations, namely that 
satisfy G a (x)W = 0. Here G a (x) = — D® c 5A t^ , is the generator of an infinitesimal gauge 
transformation, with local Lie algebra, [G a (x), G b (y)} = S(x -y)f abc G c (x), and (D^X) a = 
d tx X a + gf^A^X is the gauge-covariant derivative in the adjoint representation. 
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The proof relies upon the decomposition of -Hfp, 
H FP = H inv -(v,G)^ 



Hi nv = J 



d^X 6 



h- 



6A°(x)L 8A°{x) 6A°(x). 

= - J d 4 x = J* x (o.vrj^ (3 ' 6) 

-(«,G)t = (G,«) = jd'x J^(D^r 

where f is the adjoint with respect to the inner product defined by J dA, and (v,G) is 
the generator of the local gauge transformation v a (x). Note that H inv is a gauge-invariant 
operator, [G a (x), Hi nv ] = 0, that has exp(— Sym) as a null vector, H- 1Ilv exp(— Sym) = 0. 
Let P(A) be the normalized solution of H-ppP = for given v. It is sufficient to show that 
(W) = J dA W(A)P(A) is independent of v a (x) for gauge-invariant observables W. Let 
Sv a (x) be an arbitrary infinitesimal variation of v a (x). The corresponding change in P(A) 
satisfies SHppP + Hpp5P = 0, where SHpp = (G, Sv), so 

SP = - 5H FP P. (3.7) 

Note that SHpp P has the form of a divergence, so it is orthogonal to the null space of Hpp. 
This change in P induces the change in expectation value 

6{W) = J dA5PW 

= - I dA (tfpp 1 5H F p P) W 

r 1 t t (3 - 8) 

= - J dAP[SHl P (Hip)- 1 W] 



J dAP[(dv,G) (Hip)' 1 W], 



+ + 

where H FP = Hj nv — (v, G). It is sufficient to show that 5(W) = 0. The proof is almost 
immediate, but we must verify that the dependence of v a (x;A) and Sv a (x;A) on A does 
not cause any problem. Recall that W is gauge invariant, G a (x)W = 0, so we have 
HppW = H^ nv W, which implies that Hp P W is gauge invariant, 

G a (x) H\pW = G a (x) Hl v W = 0. 
13 



It follows by induction that (Hp P ) n W = (Hj nv ) n W is gauge invariant for any integer n, 
G a (x) (Hl P ) n W = 0, which implies that for any analytic function, /(H^W = f(H^ Y )W 
is gauge invariant G a (x) f(Hp P ) W = 0. This holds in particular for f(z) = 1/z, and we 
have G a (x) (^}p) _1 W = 0. This implies that S(W) = as asserted. Note also that if 
v = 0, then the formal solution is P = iVexp(— Sym)/^- 

This proof does not rely on Faddeev-Popov gauge-fixing which would require a gauge 
choice that selects a single representative on each gauge orbit. The Gribov critique is by- 
passed, and Singer's theorem [29] does not apply. Gauge equivalence is a weaker condition 
than gauge fixing, but sufficient for physics. In the present approach we do not attempt to 
eliminate "unphysical" variables and keep only "physical" degrees of freedom. Rather we 
work in the full A-space, keeping all variables, but taming the gauge degrees of freedom by 
exploiting the freedom of gauge equivalence. It is the unphysical degrees of freedom that 
provide a long-range correlator, and a strong candidate for a confining potential. 

Another way to obtain a gauge-equivalent probability distribution is by gauge trans- 
formation. If our class of gauge-equivalent probability distributions P v (A) is large enough, 
then it is possible to absorb an infinitesimal gauge transformation <L4 M = D M e by an ap- 
propriate change 5v of v, P V (A + D^e) = P v+ s v (A). This is true and leads to a useful 
Ward identity that is derived in Appendix C. 

There remains to choose v so it has a globally restoring property. An optimal way to 
do this is to require that the force D^v, that is tangent to the gauge orbit, points along the 
direction of steepest descent, restricted to gauge-orbit directions, of a conveniently chosen 
functional. For the minimizing functional, we take the Hilbert norm-square, F(A) = 
\\A\\ 2 = J d 4 x \A\ 2 , and we consider a variation 5A^ = r/D^v that is tangent to the gauge 
orbit in the ^-direction, where r\ is an infinitesimal parameter. We have 

5F = 2(A, 5 A) = 2 V (A, Dv) =2rj [ d A x A® {d^v a + / a6c A^ c ) 

, J ' (3-9) 

= -2rj d 4 x d^Al v a . 

Thus the direction of steepest descent, restricted to gauge orbit directions, is given by 

v a = a" 1 ^, (3.10) 

where a > is a positive constant. We shall take this optimal choice for v, so the total 
drift force that appears in the diffusion equation is given by 8 

K a x (x; A) = D™F* X + a^D^d^. (3.11) 



An alternative choice suitable for the Higgs phase was proposed in [50]. 
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Here a is a dimensionless gauge parameter. This completes the specification of the time- 
independent stochastic quantization. 

The drift force ar 1 D\d IJ ,A lJL tends to concentrate the probability distribution P(A) 
close to its region of stable equilibrium, especially if a is small. Let us find the region of 
stable equilibrium. From (3.9) we see that SF < unless A satisfies d^A^ = 0. This defines 
the region of equilibrium, which may be stable or unstable. The region of (local) stable 
equilibrium is determined by the additional condition that the second variation be non- 
negative 5 2 F > 0, for all variations 5A tangent to the gauge orbit, namely 5A = D^e, for 
arbitrary e a (x). We have just found that the first variation is given by 5F = — 2(e, d^A^). 
So we have, for the second variation, 5 2 F = — 2(e, <9 M (L4 M ) = — 2(e, d^D^e). Thus the region 
of stable equilibrium is determined by the two conditions <9 M A M = and — d fI D fJ- (A) > 0, 
namely transverse configurations A, for which the Faddeev- Popov operator — d^D^{A) is 
positive. These two conditions define the Gribov region. We expect that in the limit a — > 0, 
both conditions will be satisfied. This is the Landau gauge, with probability restricted to 
the interior of the first Gribov horizon. 

So far we have discussed Euclidean quantum field theory, which is characterized by el- 
liptic differential operators. However the above considerations also apply to the Minkowski 
case. Here the formal weight is Q(A) = iVexp^YNi], where Sym = (— 1/4) / d A xF^ u F^ v 
is the Minkowskian Yang-Mills action, where indices are raised and lowered by the metric 
9\n = 9 X>I = diag(l, 1, 1, —1). Expectation- values of gauge-invariant time-ordered observ- 
ables, are given by the Feynman path integral (W) = J dA W(A) Q(A), with (1) = 1. 
Instead of eq. (3.5), we take gauge-equivalent configurations that are solution of the equa- 
tion 



HmQ = 0, 



(3.12) 



where Hm is the corresponding Minkowskian "hamiltoniah 



n 



Hm = 



I d 4 x (i) - - , 
J SA K (x) 



9k\ (iti) 



SA x (x) 



5 



+ K\x-A)] 



(3.13) 



and the "drift force" is given by 



(x-A) 




(3.14) 



D U F» X + a~ 1 D x d^A 



The linear part of this force is 



d K g^{d»A u - d v A„) + a~ x d v g^d k A 
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which, for a > defines a regular hyperbolic operator that is invertible with Feynman 
boundary conditions. As above, one may show that the solutions to this equation for 
different values of the gauge-parameter a are gauge equivalent to each other, and for 
a — > oo, one regains the formal weight N exp(iSYM)/h- 

The drift force D\d- A is not conservative, so one cannot write down an exact solution 
to the time-independent Fokker-Planck equation HppP = 0. Nor can one express the 
solution as a functional integral over a local 4-dimensional action. However we shall, by 
successive changes of variable, transform this equation into an equation of Dyson- Schwinger 
type that may be used for perturbative expansion and non-perturbative solution. 

4. Quantum effective action in stochastic quantization 

The partition function Z(J), which is the generating functional of correlation functions 
with source J, is defined by 



It is the fourier transform (with respect to iJ x ) of the probability distribution P(A), and 
satisfies the fourier-transformed time-independent Fokker-Planck equation, 



Here K x (h-^j) is the local cubic polynomial in its argument h-^j that is defined in (3.11). 
We set Z(J) = exp\W(J)/h], where the "free energy" W( J) is the generating functional 
of connected correlation functions, in terms of which the time-independent Fokker-Planck 
equation reads 




(4.1) 



J x J x -K x (h—) Z{J) = 0. 



(4.2) 




(4.3) 



The quantum effective action 



r(A c i) = J x A ch 



X 



W(J) 



(4.4) 



is obtained by Legendre transformation from W(J), by inverting 



dW 



h dZ 
ZdT T 



(A x )j, 



(4.5) 



dJ x 
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to obtain J x = J X (A C \). In the following we shall write T(A) instead of T(A C {) when there 
is no ambiguity caused by using the same symbol for the quantum Euclidean field A and 
the classical source A = A c \. The gluon propagator in the presence of the source J is given 
by 

V xy (J) = h^i (A x - (A x )j) (A y - (Ay}j) )j 

__ dAy d 2 W (4-6) 
dJ x dJ x dJ y ' 

We note in passing that the gluon propagator V xy (J) in the presence of the source J is 
a positive matrix, since one has, for any f x , J2 xy fx^xy(J)fy = X 2 )j > 0, where 
X = J2 X fx(A x — (A x )j), which is positive since it is the expectation- value of a square. It 
is expressed in terms of the Legendre-transformed variables A and T(A) by 

= (4.7) 

Expectation-values of functionals O = O(A) are expressed in terms of Z(J), W(J) or T(A) 
by 

±\ z 



(0}j = z- 1 o(h 



dJJ 

= °{^^> m 

(0) A = 0(A + hV(A)-^y, 

where the subscript indicates that the expectation-value is calculated in the presence of 
the source J or A. In the last line, the argument of O is written in matrix notation, and 
reads explicitly A x + %V xy (A)-^-. 

The gluon propagator V xy (J) is a positive matrix, as is its inverse V~ 1 xy (A), so both 
W(J) and T(A) are convex functionals. Physics is regained when the source J is set to 0, 
namely J x = -^j- = 0. Since T(A) is a convex functional, the point -Jjj- = is an absolute 
minimum of V. In the absence of spontaneous symmetry breaking, this minimum is unique 
and defines the quantum vacuum. Thus physics is regained at the absolute minimum 
of T(A), which justifies the name 'quantum effective action'. 

In terms of the Legendre-transformed variables, the time-independent Fokker-Planck 
equation (4.3) reads 



dT 
dAZ 



dT +K X (A + hV(A) A) i] = o. (4.9) 
Here V(A) is expressed in terms of T(A) by (4.7), and K X ^A + hV(A)J^jl is evaluated 



next. 
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5. Quantum effective drift force 

We call 



Q x (A) = K x (A + hV(A)—y. (5.1) 



d_ 
OA 

the 'quantum effective drift force'. It is the expectation- value 

Q x (A cl ) = (K x ) Acl (5.2) 

of the drift force (3.11) in the presence of the source A c \, as one sees from (4.8). To evaluate 
it, we expand K X (A) in terms of its coefficient functions 

K X (A) = Kg Ay + (2l)- 1 K x %A y A z + (3!)" 1 K x % w A y A z A w . (5.3) 

The coefficient functions are found from (3.11), and are given in the explicit notation by 



ab 

fx. 1l) 

(5.4) 



= S ab [ (d 2 5 KX - d K d x ) + a~ l d K d x }6(x - y) 



abc , / 

~ 5 ym k a^' y> z ) = 9f a ° [d[\5(x - y)5 IM]K S(y - z) + d^S(y - z)5 k]x S(z - x) 

(5.6) 

+ d [K 5(z - x)5 x]fl 8(x-y)) 
Kg )a ^(x; y, z) = gf abc d^8(x - z)5 x]k 6(x - y) (5.7) 

(3) a6cd ( 4 ) abed 

= ~9 2 { F be f cde S^K]x + f ace f Me S K[x 5 u]lI (5.8) 
+ f ade f cbe S^5 x]u ) 8(x - y)5(x - z)8(x - w), 

where 5 K ^ X 5 V ^ = 5 kX 5v^ — S Kiy 5 Xfl etc. The contribution to each coefficient from 
Sym is symmetric in all its arguments, including the first. Thus K^ ab x cd v (x,y, z,w) is 
symmetric under permutations of its 4 arguments. Moreover K^ ab x (x,y) is manifestly 

i<2\ abc 

symmetric in its arguments. On the other hand the first argument of K^ t rX/x {x; y, z) is 
distinguished. 
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The evaluation of the quantum effective drift force, Q X (A) = K X (A + KD-^)1, is 
straightforward. By substitution into (5.3) we have 

d 



Q X {A) = Kg) Ay + (2!)" 1 K(%[A y + hV yu —)A z 



+ (3D" 1 K$ zw (A y + nV yu ^J (A, + KV ZV1 ^)A W ^ 
= KgjA y + (2\)- 1 K x f yz (AyA z + nV yz ) 



+ (3!)" 1 K x % w (A y + %V yu —yA z A w + hV zw ). 

Use of the identity, 

dV zw (A) = _ Vzs{A)Vwt{A) ^T(A) 

8A r zs{ ' wt{ ' 3A r dA s dAt v ; 

that follows from (V~ 1 ) Z)W (A) = q A T q^ , gives the formula for Q X (A) that is the first 
equation of next section. 



6. Basic equations for Q and V 

The first basic equation of the present method is the formula, just derived, for the 
quantum effective drift force, 

Q X (A) = K X (A) + » (2l)-'Jf(^O w + ft (2!)" 1 Jff,>„P„ : 4„ 

" * ^ Kfl M v M J^L., (6 - 1) 



where V = V(A) is the gluon propagator in the presence of the source A, and is expressed 

d 2 T{A) 
dA z dA % 



in terms of T(A) by (V x ) z ,w = qj^qa ■ This equation is represented graphically in fig. 1. 



2 

The terms of order Ti and ft correspond to one and two loops in the figure, and we write 

Q x = K x + hQ lloop , x (T) + h 2 Q 2looPtX (T). (6.2) 

The second basic equation of the present approach is obtained by writing the time- 
independent Fokker-Planck equation (4.9), satisfied by the quantum effective action, T, 
in terms of the quantum effective drift force, Q X (A), 



dT 

dT x 



<n +Q X (A) =0. (6.3) 



OA 
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This equation is of classic Hamilton-Jacobi type, with energy E = 0, and hamiltonian 
H(p, A) = Px \p x + Q X (A)}. 

The pair of equations (6.1) and (6.3) forms the basis of the present approach and 
allows a systematic calculation of the correlation functions. Equation (6.1) resembles the 
DS equation for the gluon field in Faddeev-Popov theory namely -J^- = -J=^(A + hV-Jj^)l, 
where S = Sym + Sgf + 5gh 5 and Sg^ is the ghost action. Indeed the same expressions 
appear in both equations, as is seen most easily from fig. 1, except that the contribution 
from the ghost action, ^^-(A + hVj^)l, is replaced by the term proportional to a -1 in 
the gluon vertex K^ 2 \ 

In the functional equations (6.1) and (6.3), satisfied by Q X (A) and T(A), A is a dummy 
variable, and each of these functional equations represents a set of equations satisfied by 
the coefficient functions that appear in the expansions in powers of A, 

Qx{A) = Q^Ay + (2!)" 1 Q x %, V2 A yi A y2 + ... (6.4) 

T(A) = (2\)-'T y % 2 A yi A y2 + (Vy'r^A^A^ + ... , (6.5) 

where is the proper n-vertex. The individual equations for the coefficient functions 
are conveniently obtained by differentiating (6.1) and (6.3) n times with respect to A z , 
and then setting A = 0. 

We now come to an important point. The time-independent Fokker-Planck equation 
(4.9) satisfied by T(A) is equivalent to the pair of coupled equations (6.1) and (6.3) that 
is satisfied by the pair T(A) and Q X (A). Indeed every solution of (6.1) and (6.3) yields 
a solution of (4-9) and conversely. This remark is the key to transforming the time- 
independent Fokker-Planck equation into an equation of DS type. For it turns out that 
the Hamilton-Jacobi equation (6.3) may be solved exactly and explicitly for the coefficient 
functions of T(A) in terms of the coefficient functions Q x m ^ of Q X (A), where m < n. In 
fact we shall obtain a simple algebraic - indeed, rational - formula for r' n ' = (Q) for 
every n. This allows us to change variable from the quantum effective action, V = T(<2), 
to the quantum effective drift force, Q x . It will be the last in our series of changes of 
variable, P(A) -> Z(J) -> W(J) -> T(A) -> Q X (A). 

Neither the Hamilton-Jacobi equation (6.3) nor its solution V = T(Q) contains h. 
When the solution of (6.3), V = T(Q), is substituted into (6.2), one obtains an equation of 
the form 

Q x = K x + nQ lloop , x (Q) + ft 2 g 2 ioo P ,*(Q). (6.6) 
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This is an equation of DS type for the quantum effective drift force Q x . By iteration, 
it provides the ft-expansion of Q x . The zeroth-order term K x: given in (3.11), is a local 
function of A. We shall find an approximate, non-perturbative solution of this equation. 
But first we must find T(Q). 



7. Solution for quantum effective action T(Q) 

In this section we solve (6.3) for the coefficient functions T^ 2 ) = r^ 2 )(Q) and T^ 3 ) = 
V( 3 \Q). The solution for 1?( 4 ) and for arbitrary n is found in Appendix D. 
The solution for T( 2 \ reads simply 

Note: By definition, T x n ^ X2 ___ Xn is symmetric in its n arguments, whereas Q x n i;x 1 J,... Xn has 
a distinguished first argument and is symmetric only in the remaining n — 1 arguments, 
so in general the equation V x n ^ X2 ,__ Xri = — Q x n 1 - x 1 2 \... Xn would not be consistent. However 
symmetries in fact constrain Q^- X2 to be symmetric, Qx}-,x 2 = Qx^xu as we will see, so 
(7.1) is in fact consistent. 

To prove (7.1), we differentiate (6.3) with respect to A Xl and A X2 , and obtain, after 
setting A = 0, 

4% (T x % + Q x %) + On <-> x 2 ) = 0. (7.2) 

To solve this equation for 1?( 2 ), we diagonalize all the matrices by taking fourier transforms. 
In the extended notation this equation reads, 



(7.3) 



+ [(x 1 , ax) <r* (x 2 , /U 2 , a 2 )] ) = 0, 

and we take fourier transforms, 

Q {1) t^v) = W J d A keMik-(x-y)} Q^k) 
r (2) t(^) = 6ab Wj d * ke M*k ■ (x - y)} f <J(fc). 



(7.4) 
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Color, translation, and Lorentz invariance, and use of the transverse and longitudinal 
projectors, P\^(k) = (5\^ — k\k^/k 2 ) and (k) = k\k^/k 2 give the decomposition of 
these quantities into their transverse and longitudinal invariant functions, 

Qi>) = Q {1)T (k 2 ) P A » + Q^(k 2 ) P A » 
fg(fc) = T(k 2 ) (fc) + a - 1 L(fc 2 ) P^k). 

The coefficient a -1 is introduced here for later convenience. In terms of the fourier trans- 
forms, (7.3) reads 

f [f!2 2 (*o + Q^lm +(fc,/ii,/* 2 ^-fc,/i 2 ,/ii) = o. (7.6) 

Color and Lorentz symmetries, as expressed in (7.5), constrain Q^(k) to be a symmetric 
tensor that is even in k, Q^(k) = Q^(k) = Q^(—k), as is f^ 2 J(/c). Products of such 
tensors have the same property, and as a result, the two terms in (7.6) are equal, and we 
have 

f £>,(*) l^lW + Q^Ak)} =o, (7.7) 



which proves the assertion (7.1). For future reference, we note 

(7.8) 



(D- 1 )x,(k)=t^(k) = -Q^(k) 



= T(k 2 )Pl(k)+a- 1 L(k 2 )P^(k). 



Here we have introduced the usual gluon propagator, with sources set to 0, D xy = 
V xy (A)\ A =o. It is given in terms of V by (D' 1 )^ = d ^g Ay U=o = ^xy- 

We next find 1?( 3 ). For this purpose we differentiate (6.3) with respect to A Xl , A X2 
and A X3 , and obtain, after setting A = 0, 



Ifi (T X % X3 + Q X % X3 ) + (cyclic) = 0, (7.9) 
where we have used T^ 2 ) = — Q^ 1 -*, and (cyclic) represents the cyclic permutations 

(2) 

of (1,2,3). A novelty of the stochastic method is now apparent. For Q x - X2yX3 , unlike 
Qx-L 2 i 1S n °t completely symmetric in all its arguments as it would be if the drift force 

(3) (2) 

were conservative. As a result, the equation rx,i 2 ,x 3 + Qx-x 2 ,x 3 = has no solution. This 
is already apparent to zero order in h, where Q x % 2 ,x 3 = K x 2 l 2jX3 , but K x 2 l 2jX3 is not 
symmetric in its 3 arguments, as noted above. 
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To solve (7.9) for r x x 2 i 3 , we again diagonalize the matrix T XlX by Fourier transfor- 
mation. To do so, we write the last equation in the extended notation, 

/ ^(^(x^x) [r( 3 );r; 3 (-^ 2 ^3)+Q (2) ;r;3(-;- 2 ^3)] + ) =o, 

(7.10) 

and take fourier transforms, 

Q^l\ a ^ 3 ( x ^ x ^ x 3) = ( 27r )~ 8 J d 4 M 4 M%exp(z/ci ■ x x + ik 2 ■ x 2 + ik 3 ■ x 3 ) ^ 

x 5(h + k 2 + k 3 ) & 2) «\%%(h; k 2 , k 3 ) 

r^Mi^s^ 1 '^ 2 '^ 3 ) = ( 27r )~ 8 / d 4 &id 4 &3d 4 &3exp(z/ci • xi + ik 2 ■ x 2 + ik 3 ■ x 3 ) 

x S(h + k 2 + k 3 ) f (fci, fca, fc 3 ), 

where Q (2) ^%%(/ci; fc 2 , fo) and f {3) %\%%(ki, k 2 , k 3 ) are defined only for fci + fc 2 + fc 3 = 0. 
This gives 

f P (3) -JK(fci, *a, fca) + Q (2) ^(*i; to, to)] + (cyclic) = 0. (7.13) 

We use the symmetry of f ^^2% * n ^ s ^ nree arguments to write this as 

[f £Ufci)<W 2 <W3 + (cyclic)] f ^^(fei, to, to) = - H^;\ii%(h, k 2 , k 3 ), (7.u) 

where 

H {Z) l\%\l\{ku to, k 3 ) = r^(h) Q {2)a ^(k i; k 2 , k 3 ) + (cyclic). (7.15) 

To complete the diagonalization of f^(fc), and solve (7.14) for f ^ 3 - ) ^ 1 1 %^ 3 3 (to, to, to), 
we apply a transverse or longitudinal projector to each of its three arguments, and use 
the transverse and longitudinal decomposition of T x ^(k) given in (7.8). One obtains 
r( 3 ^^f 3 (fci, to, k 3 ) in terms of its transverse and longitudinal projections, defined by 
= P^(k)X u (k) and Xj(fc) = P^{k)X v {k), 

r (3)m ^(fei,fe,*3) = " [T(fc?) +T(Jfe 2 2 ) +T(fc 3 2 )]- 1 ^ m ^(*i,fe,*3) 

(7.16) 

r^r^i^, to) = - [a" 1 ^ 2 ) + T(fc 2 ) +T(fc 3 2 )]- 1 H( 3 ) LTT »^( fc i»*2» fc 3) 

(7.17) 

etc. The corresponding formulas for 1?( 4 ) and r( n ) are found in Appendix D. 
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8. Dyson-Schwinger equation for the gluon propagator 

We have solved the second basic equation (6.3) for the coefficient functions T^ n \ 
and expressed them in terms of the Q^ m \ for m < n. We now turn to the first basic 
equation (6.1), and derive the equations for the coefficient functions by the same 

method of taking derivatives and setting A = 0. To derive the equation for we 
differentiate (6.1) with respect to A y , and obtain, after setting A = 0, 



LX 2 



ni 1 ) — /H 1 ) _ % f2n -1 d n r (3) + h f2n _1 n 

^x;y ~ ly x;y 11 K^-J 1 y x;x x ,x 2 x i Vi ^x 2 y 2 L yiy 2 y ^ \ > I ^xx 1 x 2 y ^xn 

— h (3!) K X X!X 2 X 3 -D Xl y! D X2 y 2 D Xaya ^ yi y 2 y 3 y (8-1) 

T II \*-) ly xx 1 x 2 x 3 ±y x 1 z 1 - L ^x 2 z 2 L z 1 z 2 z 3 ±y z 3 y 1 ±J x 3 y 2 L Vl ,i 



where we have again used (5.10). This equation is represented diagrammatically in fig. 2. 
In momentum space the coefficients (5.4) - (5.8) of the drift force read 

K {1) Z(x;y) = * ab (27T)" 4 / d'keM^k ■ (x-y)} K^(k), (8.2) 
K { % 1 ZZ(x 1 -,x 2 ,x 3 ) = r xa * az f d 4 k 1 d 4 k 3 d 4 k 3 exp(ik 1 -x 1 + ik 2 -x 2 + zk 3 -x 3 ) 

x5(k 1 + k 2 + k 3 ) i^ (2) M i^M3 (*i; fc 2, h) 

(8.3) 
(8.4) 



^Ci^Z^ 1 '^ 2 '^ 3 '^ 4 ) = ( 2yr ) 12 / d 4 k 1 d 4 k 3 d 4 k 3 d 4 k 4 exp(i^2k i ■ x^ 

J i=i 

x 5(h + k 2 + k 3 + k 4 ) K^ZX^, 

where 

- k^(k) = [ (k 2 5 Xil - k x k^) + a- 1 ^ ] (8.5) 

ViM2/^3 (^1) ^2 7 ^3) = ~~ ^YMAiiM2M3 (^1) ^2 ; k$) + a ^gt^iM2^3 (^1) ^2 7 ^3)7 

~ ^YMMiM2^ 3 ( /c l 5 /c 2,^3) = ig [ (^l)[ M 2^ 3 ]Mi + (cyclic) ] (8.6) 

~ f2l 

K y J ^ 2tJL3 {krMM) = ig[ (/c 3 ) Ai3 ^iM2 - (2 <-> 3) ]. 



_ ^(3)aia 2 a 3 a4 2 / fsi02era 3 a 4 e e e , faia 3 e ra 2 a 4 e e c 

4- / , oia 4 e <-a 3 a 2 e r r \ 

' J J "^l[^ 3 U ^2]M4 /• 



(8.7) 
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With = — D 1 , we obtain finally the DS equation for the gluon propagator 

- 5 ab (D^xM = - S ab [(k 2 5 KX - k K k x ) + a-'Kkx] 

_ hf aa ia2 ( 2 ! )-i(2 7r )- 4 | dh k^fc-kuh-k) 

xDx^ih) Dx^ik-h) f< 3 )»^(fci,fc-fci,-fc) 
+ h (2!)- 1 (27r)- 4 | dh K^l^ Dx^ih) 

(8.8) 

where the two-loop term is given by 

* ab Q^x,(k) = -h 2 (3!)- 1 (2tt)- 8 J dhdk 2 K^ZTxlt Dx^(h) D X2 , 2 (k2) 

x Dx 3 ^{k -h- k 2 ) t^l\% a ^{k u k 2 , k-h- k 2 , -k) 
+ h 2 (2!)- 1 (2 7 r)- 8 | dk,dk 2 K^Z\%% Dx^ih) 

x Dx 2V2 {k 2 ) r^iliikuk^-h-h) D^ih + h) 

x Dx^ 2 (k -h- k 2 ) r^^ih + k 2t k-k 1 - k 2 , -k). 

(8.9) 



9. Truncation scheme 

To obtain a non-perturbative solution of the DS equations, it is necessary to truncate 
them in some way. Needless to say, truncation remains an uncontrolled approximation 
until it is tested by varying the scheme, or by comparison with numerical simulation, as 
discussed in the Introduction and Conclusion. Moreover the truncation scheme is gauge 
dependent. This situation is familiar in atomic physics where bound state calculations 
are done in the Coulomb gauge. We shall ultimately solve the truncated system in the 
Landau-gauge limit. 

As a first step we neglect the two- loop contribution in eq. (8.8). We shall also not 
retain the tadpole term, which in any case gets absorbed in the renormalization. The 
3-vertex that we will obtain 

r (3) ™( fc i> **. h) = r a ^ f g^fo, **, fca), (9.1) 
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defined for k\ + k 2 + k 3 = 0, has the color dependence that allows us to use the identity 
f^ia 2 fa ia2 b = NS ab for $U (N) color group. As a result, the DS equation (8.8) simplifies 
to 

(t)- 1 )xM= (k 2 5 K x-k K kx)+a- 1 k K k x 

+ hN (2!)- 1 (2tt)- 4 J dh kfl^-kuk^-k), (9.2) 
xD Xl ^(ki) Dx 2 ^(k-h) tf^{k u k-k u -k). 

We convert this into a DS equation for the invariant propagator functions T(k 2 ) and 
L{k 2 ). The gluon propagator is given by 

To get the DS equation for T(k 2 ), we apply projectors Pj I/ (/c) to both free indices of (9.2), 
and obtain [P T (/c) J D- 1 (/c)P T (/c)] AAt = T(k 2 )P^(k) on the left hand side. We take the 
trace on Lorentz indices in d space-time dimensions, and use P xx {k 2 ) = d — 1, to obtain 
the DS equation for T(/c 2 ), 

T(k 2 ) = k 2 + 2{d ™ {2n)d J d d h [I T > TT {k u k) + 2I T ' TL ( kl , k) + I T ' LL ( kl , k)}, (9.4) 
1 {kl > k) ~ T{k\) T(k 2 ) 1 J 

jT,LL (k m ,2 g^gAgj -fcl. ~fc2) f^J^, fc 2 , -fc) 

ki = k — k\, and the transverse and longitudinal projections are defined in sec. 7. 

Similarly, to get the DS equation for L(k 2 ), we apply projectors (k) to both free 
indices of (9.2), and obtain [P L (k)f)- 1 (k)P L (k)}x f , = a~ x L{k 2 )P^(k) on the left hand 
side. We take the trace on Lorentz indices in d space-time dimensions, and use P xx (k 2 ) = 1, 
to obtain the DS equation for L(k 2 ), 

a- 1 L(k 2 )= a ~ 1 k 2 + ^y i Jd d k 1 [I L ' TT (k 1 ,k) + 2I L ' TL (k 1 ,k) + I L ' LL (k 1 ,k)}, (9.8) 
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where 



where 



K (2) \J \ (k;-k u -k 2 ) f^J l Wjfc fc m 





T(fc?) T(fc 2 ) 








-k) 


a — 


T(fc?) L(fc 2 2 ) 




. 2 ^U^; 




-k) 



I L ' LL {k^k) = a 2 AAlAzV ' ''^rnJ 1 2 • t 9 ' 11 ) 

The vertex iO 2 ) is given in (8.6). To complete the truncation scheme and obtain 
closed equations for T(k 2 ) and L{k 2 ), we need an approximation for the vertex T^ 3 \ We 
will approximate f ( 3 ) by its value to zero-order in fi. This vertex is expressed linearly in 
terms of by the exact formulas of sec. 7, which may be written = MQ( 2 \ where 
M = M(D). At tree level, is given by 

Q(2) = £<2) 

-_ 50) +a - 1 k (2) (9 ' 12) 

where we have used (6.2) and (8.6). Each of these terms contributes additively to f ( 3 ) = 
MQ( 2 \ Moreover Syy, is symmetric in all its arguments. As a result, it contributes 
unchanged to f( 3 ), as one sees from (7.9), and we have 

HS^Ms^i'^,^) = S^l Ll ^ 1 ^ 3 (k 1 ,k 2 ,k 3 ) + f2 ) AtlAt2M3 (/ci,/c2,/c3), (9-13) 
where fg } = MK^ } is obtained from (7.15) - (7.17) by the substitutions 

££U 3 h) - a- 1 ^?^^ (h; k 2 , k 3 ) 

rifi^a^i'^,^) -> f^Vw {h,k 2 ,k 3 ). 

~ (3 ) 

Finally, to obtain T y gt ' to zero-order in h, we substitute the tree-level propagators 

T(k 2 ) -> /c 2 , L(/c 2 ) -> fc 2 , (9.15) 

into the formulas of sec. 7. This is done in Appendix E, and gives for the vertex f^j/, 
^lll(^k 2 ,k 3 )=0 

f^lli 3 (kuk 2l k 3 ) = -ig ak2 k l~ k f +k2 (**)«, [P T (fci)P T (MW 2 
^l^fa^M = - ^9 ( ak f~ k flk 2 i pT ( k ^ ph ( k ^U, 3 " (2 - 3) ) 

rg^JiaCfci,**, **) = - ^ ( j^^q (*0/*i [^2)^3) W 3 + (cyclic)), 

(9.16) 
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(9.14) 



valid to zero-order in %. Because of the denominators, the vertex f 1 ^/ is non-local even to 
this order. Equations (9.13) and (9.16) complete the specification of that appears in 
the truncated DS equations (9.4) and (9.8) for the 2 invariant propagator functions T(k 2 ) 
and L(k 2 ). 

In Faddeev-Popov theory there are, by contrast, 3 invariant propagator functions, 
namely, these 2 plus the ghost propagator. However in Faddeev-Popov theory, the Slavnov- 
Taylor identity in its BRST version implies that the gluon self-energy is transverse, so there 
are finally only 2 independent invariant propagator functions in Faddeev-Popov theory also, 
namely, the transverse part of the inverse gluon propagator and the ghost propagator. 9 
In the present theory, the longitudinal part of the gluon propagator replaces the ghost 
propagator as the second invariant propagator function. There is no BRST symmetry in 
the present theory, but it possesses a Ward identity, derived in Appendix C, that expresses 
the effect of a gauge transformation and constrains the form of divergences. 



10. Landau gauge limit 

We now specialize to the Landau gauge limit a — > 0. We cannot directly set a = in 
the DS equations (9.4) and (9.8) because both vertices contain terms of order a -1 . With 
the gluon propagator given by (9.3), we take as an Ansatz that the invariant propagator 
functions T(k 2 ) and L(k 2 ) remain finite in the limit a — > 0. This accords with the behavior 
obtained in [6] by a Born-Oppenheimer type calculation. At a = 0, the propagator is indeed 
transverse, which is the defining condition for the Landau gauge, and L(k 2 ) does drop out 
of the propagator. However the vertices contain terms of order a -1 , and, remarkably, 
the longitudinal propagator function L(k 2 ) does not decouple at a = 0, but remains an 
essential component of the dynamics! 

We next determine the a-dependence of the vertices asymptotically, at small a. By 
(8.6), we have = — Syy[ + a~ 1 K^\ so this vertex contains a term of order a and 
a term of order a -1 . We take the asymptotic limit of (9.16) at small a, and obtain the 
interesting a-dependence 

^lll(h,k 2 ,k 3 )=o 

f g^Si^fci, fc 2, h) = a" 1 nfrX^iK k 2 , k 3 ), 

9 In practice the truncated DS equations in Faddeev-Popov theory violate the Slavnov- Taylor 
identities to some extent. 
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where 

ntllisfahM = -ig (h), 3 [P T (h)P T (k 2 )U, 2 

im 3 (k u k 2l k 3 ) = -ig (k 2 ), 2 [P T (k 1 )P L (k 3 )}, 1 , 3 - (2<-> 3) (10.2) 

K 2 K 3 

t^ s {kiMM) = -ig [^2)^3) Wa + (cyclic) 

are independent of a. These quantities are anti-symmetric in their three arguments so, for 
example, 7^ JJ^fci, k 2 , fa) = -7J J 3 J 2 (fci, fc 3 , fc 2 ), etc. We see that f( 3 ) also contains a 
term of order a and a term of order a -1 . 

The DS equation for L(k 2 ), eq. (9.8), is consistent with our Ansatz in the Landau 
gauge limit only if the leading term on the right is also of order a -1 . This is non-trivial, 
because both vertices contain terms of order a -1 , so in principle terms of order a -2 could 
appear on the right hand side which would invalidate our Ansatz. 

We now derive the DS equation for L(k 2 ) in the Landau-gauge limit by evaluating in 
succession the terms (i) J L ' TT , (ii) J L ' LL ? and (iii) I L > TL that appear on the right hand 
side of (9.8), in the limit a — > 0. 

(i) Consider eq. (9.9) for J L ' TT . It contains no explicit powers of a. Moreover the 
vertex r^^T^ k 2 , — k), given in (10.1), is of order a . Thus the inconsistency of a 
term of order a~ 2 is avoided, and this intermediate state will give a contribution of required 
order a -1 only if the vertex, K^Yx^-2 ~~ ~k 2 ), gives a contribution of order a -1 . The 
term a~ 1 K^\x 1 \ 2 (k; —ki, —k 2 ) in (8.6) is in fact of this order. The projected components 
of Kg/ are easily read off (8.6) by writing = [P T (k) + P L (k)]x fI , which gives 

kf^ M {k 1 -MM = w ([P T (h)+p L (h) }[p T (k 2 )+p L (k 2 ) } ) 

- (2<->3). 

(10.3) 

The polarization vector (ks)^, 3 is purely longitudinal, as is {k 2 )^ 2 , and this implies 

K$\ZZ^ -ku -* 2 ) = i^L^; -fei, -* 2 ) = 0. (io.4) 

Thus there is no contribution of the required order a -1 from 7 L ' TT . 

(ii) Consider eq. (9.11) for / L > LL . It has the coefficient a 2 , so there is no contribution 
of the required order a -1 from J L > LL either. 
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(iii) Now consider eq. (9.10) for 7 L ' TL . It has the coefficient a. So when each vertex 
is of order a -1 there is an overall contribution to a~ 1 L(k 2 ) of the required order a -1 . As 
a result, the DS equation (9.8) for L(k 2 ) simplifies in Landau gauge to 

L(k 2 ) = k 2 + hN (27T)" 4 J d*h ^SM^tl ^^y^f 2 ^ 1 - ~*° ■ ( 10 - 5 ) 



By (10.3) and (10.2), we have 
^aTa'^ ~ku -**) llUikn fc, -k) = - 0^ [k . P T {kl) . k + k . p T {ki) . ^ 

n 2 k 2 

= -2-^-% I k-P T (k 1 )-k, 



(10.6) 

where we have used k 2 = k — k\. Note that a factor of the external momentum k appears 
at each vertex. This corresponds to the factorization of external ghost momentum in the 
Landau gauge in Faddeev-Popov theory. This gives the truncated DS equation for L(k 2 ) 
in Landau gauge, 



, 2hg 2 N f, 4? k 2 [k 2 k 2 -(k-h) 2 ] 



L{k ] ~ k J dh k 2 (k 2 + k 2 )T(k 2 ) L(k 2 )- (1 °- 7) 

The DS equation for T(k 2 ), eq. (9.4), is consistent with our Ansatz in the Landau 
gauge limit only if the leading term on the right is also of order a . This is non-trivial, 
because both vertices contain terms of order a -1 , so in principle terms of order a -1 and 
a~ 2 could appear on the right hand side which would invalidate our Ansatz. 

We now derive the DS equation for T(k 2 ) in the Landau-gauge limit by evaluating in 
succession the terms (i) J T ' TT j (ii) / T ' TL , and (iii) / T ' LL that appear on the right hand 
side of (9.4), in the limit a — > 0. 

(i) Consider eq. (9.5) for / T ' TT . It contains no explicit powers of a. By (10.1) and 
(10.4), the vertices from K^JJ^^k; —k\, —k 2 ) and f gt ^(ki, k 2 , —k) vanish, and we 
obtain from (8.6) and (9.13), 

K {2) l^ T Ak--k^-k 2 ) = - S^lZlik, -h,-k 2 ) (10.8) 

and 

f (3) J J 2 h, -k) = S&lllfa, k 2 , -k). (10.9) 
This gives a contribution of the required order. 
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(ii) Next consider eq. (9.6) for 7 T ' TL . It has coefficient a. By (9.13) and (10.1), we 

have 

r^JJfo, fa, -k) = s< 3 > ^(h, to, -k) + 7 (3) ;^J(to, fa, -ft), (10.10) 

which is of order a . Thus only that part of the vertex K^J^^^k; —fa, —fa) that 
is of order a -1 will contribute to the desired order a . However from (10.2) for 
7 ( 3) (ki, fa, fa), and by evaluation of S^^^^fa, fa, -k), one obtains 

f (3) J^J(fci,to,-fc)=o. (10.11) 

(iii) Finally consider eq. (9.7) for J T ' LL . It has coefficient a 2 . To get a net contribution 
of order a , we make the substitutions of the relevant projected vertices, 

K^lUAk; ~ki, ~fa) - a^k^X (k; -fa, -fa) 
f (S) ^Z(fa, fa, -k) - a-^\^{fa, fa, -k), 

by (8.6), (9.13), and (10.1). Again the conclusion is consistent with our Ansatz. 

We have now found all the terms on the right hand side of (9.4) that contribute to 
T(k 2 ) in the Landau gauge limit, namely, 

_ c(3)TTT/i _u _u\ 5(3) T T T(u jl _M 

I T ' TL (fa,k) =0 (10.14) 

jT,LL ( 1, , n _ j^j^Aa^j ~ fc 2) 7 (3) AiA 2 I( fc l» fc 2, ~fc) ^ ^ 

and fa = k — fa. The last term is given explicitly by 
I T > LL (fa,k) = - g 



T LL l _ 2 (fc 2 + fcf ) fa ■ P T (k) -fa-klfa- P T (k) -fa-klfa- P T (k) ■ fa 



(kf + ki) L(k{) L(kj) 

(10.16) 

where we have used fa ■ P T (k) ■ fa = fa ■ P T (k) ■ fa = —fa ■ P T (k) ■ fa. The non-local 
denominator (k 2 + k 2 )~ x has cancelled out of this expression. 

We have obtained a consistent Landau gauge limit of the truncated DS equations for 
the invariant propagator functions T(k 2 ) and L(k 2 ). As asserted, the invariant longitudinal 
propagator function L(k 2 ) does not decouple in this limit. The reader will have noticed 
a striking similarity to the corresponding equations in Faddeev-Popov theory, with the 
longitudinal propagator replacing the ghost propagator. 
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11. Infrared critical exponents 

We shall solve the the DS equations (9.4) and (10.7) for the asymptotic forms of T(k 2 ) 
and L(k 2 ) in Landau gauge at low momentum. We suppose that at asymptotically small 
k, they obey simple power laws, 

T(k 2 ) ~ C T (k 2 ) 1+aT 

(11.1) 

L(k 2 )~C L (k 2 ) 1+ «\ 

where ax and a l are infrared critical exponents whose value we wish to determine. Canon- 
ical dimension corresponds to ay = a/, = 0. As explained in sec. 3, we know that in the 
Landau gauge limit, a — > 0, the gauge field A is constrained to be transverse d ■ A = 0, 
and to lie inside the Gribov horizon, that is to say, where the Faddeev-Popov operator is 
positive, —d ■ D(A) > 0. The transversality condition is satisfied by our Ansatz. As has 
been shown many times [22], [26], the positivity condition strongly suppresses the low- 
momentum components of A(k). Recalling that the transverse part of the gluon propgator 
is given by D T (k 2 ) = {\A(k)\ 2 }, we look for a solution for which D T (k 2 ) = 1/T(k 2 ) is 
suppressed at low k, so T(k 2 ) is enhanced at small k compared to the canonical power 
T(k 2 ) = k 2 . This means a T < 0. 

We now estimate the power of k of the various terms in the DS equation (9.4) for T(k 2 ). 
The analysis is similar to the Faddeev-Popov case [2], [3], [4], [5], [6]. The left hand side 
has the power (k 2 ) 1+aT . The tree-level term is k 2 , so with ar < 0, the tree level term is 
subdominant in the infrared and may be neglected. To evaluate the loop integral f d 4 ki, 
asymptotically at low external momentum k we take k to be small compared to a QCD mass 
scale, \k\ « Aqcd, and we rescale the variable of integration according to k± = \k\x^. 
We now have a dimensionless integral in which the QCD mass scale appears only in the 
very small ratio |/c|/Aqcd- I n the asymptotic infrared limit, this ratio goes to 0, and 
everywhere in the integrand we use the asymptotic forms (11.1). This is equivalent to 
using the asymptotic forms (11.1) everywhere in the original integral. We shall see that 
the resulting integral is convergent, which means that the integral is effectively cut off at 
momentum k± ~ k. 

We now estimate the contributions of the terms J T ' TT and I T ' LL , eqs. (10.13) and 
(10.15), to the right hand side of the DS equation (9.4), by simply counting powers of k 
and k\. One finds that, after integration J d 4 ki : these terms are of order {k 2 ) 1 ~ 2aT and 
(k 2 Y~ 2aL respectively, while the left-hand side is of order (/c 2 ) 1+aT , with a T < 0. The 
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powers match on both sides only if aj, > 0. In this case, I T > LL is the dominant term on 
the right, and by equating powers of k, one obtains 

a T = -2a L , (11.2) 

and cxl > 0. We retain only the dominant term J T > LL on the right in (9.4), which simplifies, 
for arbitrary space-time dimension d, to 

Gr (k 2 )^ " ^ fd d k ^ ^ ~ ih " ^ (11 3) 

Ct {k > - (d - 1) c 2 (2^ k 2 J d kl (jfe?)i+^ (k 2 y+^ ' [iL6) 

where k<i = k — k\, and we have used (10.17). This agrees with eq. (6.14) of [6] in Faddeev- 
Popov theory. We write it as 



CtC 



(1L4) 

where It is evaluated in Appendix F. We have generalized to arbitrary space-time dimenion 
d, and we take d in the range 2 < d < 4. By equating powers of k for arbitrary d, we find 
that the critical exponents are related by 

a T + 2a L = - (4-d)/2. (11.5) 

The last integral is ultraviolet convergent provided that oll > (d — 2)/4, which corresponds 
to ar < —1- For d = 4, we obtain ct^ > 1/2 as the condition for convergence of the integral. 
Now consider the DS equation (10.7) for L(k 2 ) in the infrared asymptotic limit, 

Cr(k*V^-e 2hg2N fd*k, ^kl-{k- kl f 

CL{k ] ~ k C T C L (27T)4 J d kl (Jfe2)2+a T {k 2 )aL (^2 + ^2)' l 11 ^ 

for d = 4. By power counting, the integral on the right has the power (k 2 ) 1 ~ aT ~°' L . This 
agrees with the power on the left, provided cut = 2chl, which is identical to the previous 
equation. However we have also previously found oll > 0. In this case, the tree level term 
k 2 is dominant in the infrared, and the equation appears inconsistent. However the degree 
of divergence of the integral is 2a l, so the integral diverges for > 0, and a subtraction 
is required. The integral contains an explicit factor of k 2 , and the divergence is of the 
form Bk 2 , where B is an infinite constant. We subtract the integrand at k = 0, which 
makes the integral vanish more rapidly than /c 2 , and add bk 2 on the right, where b is an 
arbitrary finite constant. The dominant terms are now the tree level term k 2 and bk 2 . 
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For the equation to be consistent, the subtraction term must precisely cancel the tree-level 
term, so b = — 1. This gives 



C T C L (2tt) 4 7 " '" A (fc?) 2 +^ V (fc?)«* fef (fcl + A; 2 )- 

(11-7) 

This integral is also convergent in the infrared for «t = — 2a^ < 0. The right hand side now 
vanishes more rapidly than k 2 . This conclusion, agrees with the "horizon condition" [25], 
and with the confinement criterion of Kugo and Ojima in the BRST framework [27], [28]. 
Conversely we could have imposed the horizon condition on the DS equation for L(k 2 ), 
and derived the suppression of the transverse propagator 1/T(k 2 ) at low momentum. 

The subtracted expression on the right is most simply evaluated by continuing in 
space-time dimension d. In this case one can ignore the subtraction term, and evaluate the 
unsubtracted integral with dimensional regular ization for d < 4, and continue the resulting 
expression to d = 4, 

c,(k 2 )^ gViv r d k 2 k 2 -(k- kl ) 2 

CL[k) ~ c T c L {2*y] dkl (k 2 ) 2 +^ { k 2 )^ (k 2 + k 2 y UL8j 

The denominator k\ + k 2 results from the non-local expression for the vertex. One obtains 
the corresponding equation for the ghost propagator in Faddeev-Popov theory, eq. (6.15) 
of [6], from this equation by the substitution k ^ k2 — ► p-- 

By equating powers of k for general space-time dimension d, one again gets (11.5), 
and we see that the DS equations for the transverse and longitudinal parts are consistent. 
The degree of divergence of this integral is 2a l, and after one subtraction its degree of 
divergence is 2aL — 2, so the subtracted integral is convergent provided that cxl < 1, or 
equivalently that «t > — 2 — (4 — d)/2. From this and our previous bound, we conclude 
that for d = 4, this subtracted integral and (11.3) are both finite provided that «i is in 
the range l/2<a L <l, or equivalently that ar is in the range — 2 < a T < — 1. We write 
the preceding equation as 



CtC 



Tig 2 N 

where II is evaluated in Appendix F. 

Upon comparison with (11.4), we obtain 



= Il, (H.9) 



I T (a L ) = I L (a L ) (11.10) 
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which determines the critical exponent a^. From Appendix F, this gives, for d = 4, 

r 2 (2 - a L ) T(2a L - 1) _ 3 + 2a L + 2) T 2 (l - q L ) r(2a L + 1) 

r 2 (l + a L ) T(4-2a L ) a L T(a L + 2)r(a L + 3)r(2 - 2a L ) ' 1 ' ' 

Both expressions are finite and positive, as they should be, for a L in the interval 1/2 < 
«i < 1. Moreover at ct_L = 1/2, the left-hand side diverges whereas the right is finite. 
On the other hand at «l = 1, the left-hand side is finite whereas the right diverges. 
Consequently there is at least one root in the interval 1/2 < a L < 1. After cancelling 
T-functions, the last two equations give the quartic equation 

49a| - 189a| + 133a| + 117a L - 74 = 0. (11.12) 

From a numerical investigation it appears that there is only one root in the interval 1/2 < 
aL < 1, with the value 

a L w 0.5214602698 

(11.13) 

a T = -2a L « -1.04292054. 



12. Conclusion 

We derived time-independent stochastic quantization from the principle of gauge 
equivalence which states that probability distributions that give the same expectation 
values for all gauge-invariant observables are physically indistiguishable. This quantiza- 
tion is expressed by an equation for the Euclidean proabability distribution P(A) that is of 
time-independent Fokker-Planck form, with a corresponding equation for the Minkowski 
case. By making several changes of variable, we transformed this equation into an equation 
of DS type, suitable for non-perturbative calculations. The most novel of these changes of 
variable is accomplished when the equation for the quantum effective action V is exchanged 
for an equation for the quantum effective drift force Q x . We then adopted a truncation 
scheme and obtained a consistent Landau gauge limit, a — > 0, and found, remarkably, that 
the longitudinal propagator function L(k 2 ) that appears in the longitudinal part of the 
gluon propagator D L = a/L(k 2 ), does not decouple in the a — > limit, but plays a role 
similar to the ghost in Faddeev-Popov theory. 

We calculated the infrared critical exponents that characterize the asymptotic form 
at low momentum of the transverse and longitudinal components of the gluon propagator 
in Landau gauge, D T ~ l/(/c 2 ) 1+aT , and D L ~ a/(k 2 ) 1+aL , and obtained the values 
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a L « 0.5214602698 and a T = -2«l « -1.04292054. In the Landau-gauge limit a -> only 
the transverse part survives. As a function of k, it vanishes at /c = 0, albeit rather weakly, 
D T ~ (/c 2 )~ 1_aT ~ (/c 2 ) 043 . On the other hand, the longitudinal part of the propagator, 
is long range, D L ~ a/(/c 2 ) 1521 . Qualitatively similar values have been obtained recently 
for the infrared critical exponents of the gluon and ghost propagators in Landau gauge 
from the DS equation in Faddeev-Popov theory, using a variety of approximations for the 
vertex, [2], [3], [4], [5] and [6], in particular, [5] and [6], «t = = —1.1906, and 

ac = 0.595353 respectively. As we have argued recently [6], these calculations in Faddeev- 
Popov theory should be interpreted as including a cut-off at the Gribov horizon. This 
makes them similar in spirit to the present calculation for which, as shown in sec. 3, the 
probability also gets concentrated inside the Gribov horizon in the Landau gauge limit 
a — > 0. Reassuringly, the solutions of the DS equation in Faddeev-Popov theory and in 
the present time-independent stochastic method are in satisfactory agreement. 

We comment briefly on the physical significance of our results, (i) We have avoided 
gauge fixing and instead derived the equation of time-independent stochastic quantization 
from the principle of gauge equivalence, thereby overcoming the Gribov critique. Since 
we do not gauge fix, we do not brutally eliminate "unphysical" variables and keep only 
"physical" degrees of freedom, which would violate Singer's theorem [29]. Instead, we 
gently tame the gauge degrees of freedom by exploiting the principle of gauge equivalence, 
(ii) We derived a set of equations of DS type that was solved approximately but non- 
perturbatively in Landau gauge asymptotically at low momentum, (iii) The values we 
obtained for the infra-red critical exponents of the gluon propagator in Landau gauge are 
in satisfactory agreement with corresponding values in Faddeev-Popov theory, and also 
with numerical simulations, (iv) A striking result of this investigation is that the invariant 
longitudinal propagator function L(k 2 ) does not decouple in stochastic quantization, even 
though the longitudinal part of the gluon propagator D L = a/L(k 2 ) vanishes with the 
gauge parameter a in the Landau gauge limit a — > 0. Indeed, because some vertices are 
of order a -1 , transverse gluons exchange longitudinal gluons as virtual particles, with an 
amplitude that remains finite in the limit a — > 0. Thus, while ghosts are absent in time- 
independent stochastic quantization, they are replaced dynamically by the longitudinal 
part of the gluon propagator in the Landau gauge limit. In fact, the DS equations (9.4) 
and (10.7) for T(k 2 ) and L(k 2 ) bear a remarkable similarity to the DS equations for the 
gluon and ghost propagators D and G in Landau gauge in Faddeev-Popov theory, with the 
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correspondences D <-> 1/T and G <-> 1/L. In both cases, it is the ghost loop or longitudinal- 
gluon loop that gives the dominant contribution to the transverse-gluon inverse propagator 
in the infrared region, and causes suppression of the would-be physical, transverse gluon 
propagator at k = 0, a signal that the gluon has left the physical spectrum. 

To conclude, we mention some challenging open problems, (i) The possibility of com- 
parison with numerical simulations is an essential and promising aspect of the present 
situation. Any DS calculation involves a truncation which remains an uncontrolled ap- 
proximation, without further investigation. It may be controlled by varying the vertex 
function [5], or by extending the calculation to include the vertex self-consistently. For- 
tunately, comparison with numerical simulation provides an independent control. In this 
regard we note that the stochastic quantization used here may be and, in fact, has been 
effected on the lattice in numerical simulations by Nakamura and collaborators [51], [52], 
[53], [54], [55]. A direct comparison with this data would require a solution of the DS 
equations for finite gauge parameter a, or extrapolation of lattice data to a = 0. Natu- 
rally, a comparison of numerical results with asymptotic infrared calculations also requires 
control of finite-volume lattice artifacts, (ii) Conversely the results of the DS calculations 
suggest new numerical calculations. In particular our prediction that, for small values of 
the gauge parameter a, the longitudinal part of the gluon propagator is long range, should 
be tested numerically, (iii) The present scheme is not based on a local action, but rather on 
the DS equations of time-independent stochastic quantization. Renormalizability follows 
from the indirect argument that correlators of the 4-dimensional time-independent formu- 
lation used here coincide with the equal-time correlators of a local, 5-dimensional theory 
whose renormalizability has been established [49]. The renormalization constants of the 
5-dimensional theory were calculated some time ago at the one-loop level, and were found 
to yield the usual (3- function [57] , [58] , [46] . However a direct proof of renormalizability in 
the present time-independent formulation remains a challenge. The Ward identity, derived 
in Appendix C, is a first step, (iv) One should extend the solution obtained here for the 
asymptotic infrared region to finite momentum k. (v) The Landau gauge is a singular 
limit, a — > 0, of the DS equations for finite gauge parameter a. It would be valuable to 
also solve the DS equations for arbitrary, finite a. (vi) One should extend the solution 
of the DS equations to include quarks, (vii) As we have explained, our results are in- 
tuitively transparent and lend themselves to a simple confinement scenario in which the 
would-be-physical transverse gluon leaves the physical spectrum. However it is clear that 
our discussion of confinement remains at the level of a scenario because we have dealt 
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here only with the gluon propagator which is a gauge-dependent quantity. This is only 
a first step in a program, some of whose elements have just been indicated. Clearly the 
goal is to calculate gauge-invariant quantities. Gauge-invariant states, the hadrons, ap- 
pear as intermediate states in gluon-gluon and quark-anti-quark scattering amplitudes. 
One must extend to this sector the solution of the DS equations obtained here, and of the 
Bethe-Salpeter equations that follow from them. 
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Appendix A. Time-independent Fokker-Planck equation for quarks 

We extend to quarks the derivation of time-independent stochastic quantization from 
the principle of gauge equivalence. Following the method of sec. 3, we seek a weight 
function P = P(A, ip, V>) that depends on the gluon and quark and anti-quark fields. We 
wish to establish a class of gauge-equivalent normalized distributions that includes the 
formal gauge-invariant weight P = N exp(-S) as a limiting case. Here 

S = S YM + J d 4 x i;(m+ (A.l) 

is the Euclidean action of gluons and quarks, = l^D^ = ^^{d^+gA a ^t a ), where {7^, 7^} = 
28^, and the t a are the quark representation of the Lie algebra of the structure group, 
[t a ,t b ] = f abc t c . We take P to be the solution of H FP P = 0, where we now specify the 
extended Fokker-Planck hamiltonian. 

As in sec. 3, we take H FP to be of the form 

H FP =H inv -(v,G)^ 

- H - + I d4 < - 4^ + 9 h m ) * e ' (A ' 2) 

where H inv is a gauge-invariant operator, specified below, that has exp(— S) as null vector, 
H inv exp(— S) = 0, and the Grassmannian deriveratives are left derivatives. Here 

G » = " " 9(i>) ^M + 9m JW) (A ' 3) 
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is the generator of local gauge transformations that satisfies 

[G a (x),A b x (y)] = -D a J> 8(x-y) 
[G a {x)^{y)) = - g t^(x) 8{x - y) 

[G a {x)^{y)]=g^{x)t a 5{x-y) (A.4) 
[G a (x),G b (y)] = 5(x-y) gf abc G c (x) 
[G a (x),H inv ]=0. 

With gauge-invariant observables defined by the condition G a (x)W = 0, the proof of 
sec. 3, that the expectation-value of gauge invariant observables (W) = J dAd^d^ W P is 
independent of v, applies here as well. As explained in sec. 3, we take v a (x) = a~ 1 d\A\(x), 
where a is a gauge parameter. 

There remains to specify H iriv . We suppose that it is a sum of gluon and quark and 
anti-quark hamiltonians, 

H inv = H 1 + H 2 + H 3 (A.5) 
where Hi is the gauge-invariant gluon hamiltonian as in sec. 3, 




(A.6) 



For the quark and anti-quark hamiltonians we take 




(A.7) 



where N 2 and A3 are gauge-covariant kernels with engineering dimensions of mass. All 
terms in Hpp contain a derivative on the left, which assures that Hpp has a null eigenvalue, 
for we have J dAdt^dtfj Hp-pF = for any F. The corresponding right eigenvector P that 
satisfies HppP = 0, is the physical distribution that we seek, that depends on the gauge 
parameters. Each of the operators Hi satisfies Hiex.p(—S) = 0, so also H inv exp(— S) = 0. 
This assures the applicability of the proof of sec. 3, namely that the normalized solutions 
for different v are gauge equivalent, P V (A) ~ P V /(A), and include Aexp(— S) as a limiting 
distribution. 

We have obtained this result without any assumptions about the kernels N 2 and AT 3 , 
apart from gauge covariance (and regularity). This would not be consistent unless the 
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normalized solutions for different choices of the kernels give gauge-equivalent distributions 
Pn (A) ~ Pn' (A) or, in other words, the parameters that specify N2 are and N3 are gauge 
parameters. We prove this directly. 

Let 5N2 be an infinitesimal variation of N2. It induces a corresponding change in H2 

sB *=j**h 6N *& s £)- (A - 8) 

The corresponding change in P satisfies 5H 2 P + H FP 5P = 0, so 5P = — Hpp5H 2 P. Let 
W be a gauge-invariant observable. We have 



8(W) = J dAdipdip 5P W 

= - J dAdi/jdifj H~^5H 2 P W 

= - J dAdipdip P 8h\ (Hip)' 1 W (A.9) 
= - J dAdipdip P 5h\ (h}^)- 1 W 
= - J dAdtfjdifj 5H 2 P (hIx 1 W, 

where we have used (ifp P ) _1 W = (H^) -1 W which holds for a gauge-invariant observ- 
able, as was shown in sec. 3. Moreover 5h\ (H^) -1 W is gauge invariant, so the last 
expression is independent of the gauge parameter a, as was also shown in sec. 3, and we 
may evaluate it for a — > 00. We have 

lim 5H 2 P = [ d^x^- 5N 2 f-t + lim P = 0, (A.10) 

a^oo J dip ydip dip/ a^oo 



because lim^oo P ~ N ex.p(—S). Thus 6{W) vanishes, as asserted. 
The quark action satisfies 

S qu = J d 4 x ij(m+ fl)ip = - J d 4 x ipC-\m+ p)C$, (A.ll) 

where C is a numerical matrix that acts on spinor and group indices and satisfies C _1 7 M C = 
~7* r and C~H a C = -(t a ) tr , so we have 

6 S 

-j— = (m+ 0)<p(x) 

SS (A.12) 
— = -C-\m + P)CHz), 
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where the Grassmannian derivatives are left derivatives. The most general expressions for 
N 2 and N 3 that are local, gauge-covariant, and have dimension of mass are 

N 2 = m 2 - b 2 P 

(A.13) 

N 3 = - C- x {m z -b z ip)C, 

where the and bi are gauge parameters. We expect that the kernel ultimately appears 
in the denominator in loop integrals so, to improve convergence, we should take b 2 ^ and 
63 7^ 0. For the gauge choice to respect charge-conjugation invariance, we take b 2 = 63, 
and mi = m 2 = cm, which gives 



H 2 = Jd^(cm-bP)(± + (m+ m) 

d 4 x^ {-l)C-\cm-bP)C (— + (-l)C- 1 (m+^)C^). 



(A.14) 



where b and c are gauge parameters. This gauge choice also respects chiral symmetry in the 
limit m — > 0. One may show that the eigenvalues of H 2 and H3 are the eigenvalues of fermi 
oscillators, with frequencies A n that are the eigenvalues of the operator (c m—b p)( m+ p), 
which for b = c > 0, simplifies to b(m 2 — p 2 )- In this case H 2 and H3 have the unique 
null eigenvector exp(— S), and all other their eigenvalues are strictly positive, as occurs for 
Hi. Indeed Hi satisfies 

exp( 5 /2)iW- S /2) = /^ (^-(V^X-^-d/^), (A.15) 

where the operator on the right is manifestly positive, with the unique null-vector 
exp(— S/2). Thus Hi has the unique null- vector exp(— S) : and all its other eigenvalues 
are strictly positive. However we expect that b and c must be kept as independent con- 
stants when needed as renormalization counter-terms. 

Altogether, the total Fokker-Planck hamiltonian, including quarks, is given by 



H 



FP 



J I 8A^ V 8 8 An J 

$ , , t*n ( <5 , 8S \ 



where a > 0, b > and c > are gauge parameters. 



(A.16) 
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Appendix B. Time-independent stochastic quantization on the lattice 

We briefly outline how to extend time-independent stochastic quantization to lattice 
gauge theory. To each link (x, /j,) of the lattice is associated a variable U Xyfl G SU(N). These 
variables are subject to the local gauge transformation U Xifl — > 9 U X ^ = g~ 1 U XjfI g x+ {x, 
where g x G SU(N) is associated to the site x of the lattice. Observables W(U) are invariant 
under this transformation, W( 9 U XjfJ ,) = W(U Xifl ). Expectation values are calculated by 
(W) — J dU W(U) Pw{U), where dU is the product of Haar measure over all link variables 
of the lattice, and Pw = iVexp(— Sw) is the normalized probability distribution associated 
to the gauge-invariant Wilson action Sw- 

We shall exhibit a Fokker-Planck hamiltonian H FP for the lattice, such that the posi- 
tive normalized solutions P to H FP P = are gauge equivalent to Pw, P ~ Pw- Let J£ be 
the Lie differential operator associated to the group variable on the link (x, /J,), that satisfies 
the Lie algebra commutation relations [J Xjfl , Jy^] = S xy 5^f abc J x ^- And let G x be the 
generator of local gauge transformations that is defined by (1 + J2 X £xG x )F(U) = F{ 9 U), 
where g x = 1 + e x is an infinitesimal local gauge transformation. These generators 
satisfie the Lie algebra commutation relations of the local gauge group of the lattice 
[G x ,G b y ] = 5 xy f abc G%, and may be expressed as a linear combination of the J's. A 
hamiltonian with the desired properties is given by 

H FP = H- inv -(v,G)^ 
(v,G) = ^2v x G x , 

X 

where f is the adjoint with respect to the inner product dU, and v x is a site variable with 
values in the Lie algebra. Indeed, the argument of sec. 3 holds here, with the substitution 
<Sym(^4) — ► Sw(U), and shows that the probability distributions P v for different v, defined 
by HppP = 0, are gauge equivalent to each other P v ~ P v > and to Pw- As in sec. 3, we 
choose v x (U) so the infinitesimal gauge transformation g x = 1 + e t a v x is the direction of 
steepest descent in gauge orbit directions of a minimizing functional F(U). A convenient 
choice is F(U) = J2( xy ) tr (-^ — U( xy )), where the sum extends over all links (xy) of the 
lattice. 
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Appendix C. Ward Identity 

In sec. 3 we showed that probability distributions P V (A) for different v are gauge 
equivalent, P V (A) ~ P V '(A). Another way to make gauge equivalent distributions is by 
making a local gauge transformation, because, for all gauge-invariant observables W(A), 
this cannot change the expectation value J dA W{A)P V (A) = J dA W(A)P V ( 9 A), and 
we have P V (A) ~ P V ( 9 A). If the class of gauge-equivalent distributions P V (A) that was 
introduced in sec. 3 is large enough, then the gauge transformation corresponds to a change 
of u, 

P V (°A)=P V ,(A) (C.l) 

for some v'. This is in fact the case, and provides a Ward identity. 

To prove this, we apply the infinitesimal local gauge transformation 1 + (e, G), where 
(e, G) = J d 4 x e a {x)G a {x), to the time-independent Fokker-Planck equation (3.5)and (3.6), 

[l + (e,G)] [# inv + (G»] P v = 0. (C.2) 

From the commutation relations 

[(e, G),H inY ] = 

[(e,G),G a (x)] = -r bc e b (x)G c (x), ^ 

we obtain 

[(e,G),(G,v)] = (G,5v), (C.4) 

where 

5v a = [(e, G),v a ] + f abc e b v c , (C.5) 

and, to first order in e, 

[^inv + (G, v + Sv)} [1 + (e, G)] P v = 0. (C.6) 

Note that while v a (x) and e a (x) are both local gauge transformations, v a (x) = v a (x. A) 
depends on A, but e a (x), by assumption, does not. By comparison with the defining 
equation for the probability distribution P v +sv, 

Py+Sv — 0, (C.7) 
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we conclude that the gauge-transformed probability distribution [1 + (e,G)] P V (A) = 
P V (A + De) coincides with P v +sv, 

P v (A + De)=P v+5v (A), (C.8) 

where 5v is given above. This states how a gauge transformation is absorbed by a change 
in v, and provides the Ward identity. 

This identity is inherited by the functionals we introduced, the quantum effective 
action Y v and the quantum effective drift force Q v , and we have 
F V (A + De) =T v+Sv (A) 

Q a v (x, A + De) = Q a v+5v (x, A) - f abc e b (x) Q c v+5v (x, A). ^ 
We now specialize to v = a~ 1 d ■ A, and obtain 

Sv a = a- x d ■ D ac (A) e c + a-7 a6c e & ■ A c 



= a _1 Z) ac (A) • d e c = a^fdV + f abc A b d M e c ] 



i < 

(CIO) 

y /X ^ c J- 

Only the derivative of e appears here because, for v = a~ 1 d ■ A, the probability distribu- 
tion P V (A) is invariant under global (x-independent) gauge transformations. We further 
specialize to a linear dependence of e on x, e a (x) = rj^x^, where the rf^ are infinitesimal 
constants. In this case we have 

Sv a = a~^r hc A\ r£. (C.ll) 

Although this breaks Lorentz invariance, it does not break translational invariance, so the 
perturbed hamiltonian defined by 

H FP + 5H FP = H inv + (G, v) + (G, Sv) 

(C.12) 

= H inv + a- x {G, d ■ A) + a- x (G, A^ x Vlt ), 
where (A^ x 7] fI ) a = f abc A^rj^, is translationally invariant, even though A^ + D^ c e c = 
^ + + / abc ^ Vv x vi has an explicit x-dependence. Moreover the inhomogeneous term 
a~ 1 d 2 e in 5v vanishes with this choice of e, so A = remains the classical vacuum. 
Without further calculation we conclude that the transformed quantum effective action 
T V (A + De) = T v+ Sv(A) is a translationally invariant functional of A for v = a~ 1 d ■ A and 
e a (x) = rj^x^, with A = as classical vacuum. 

More generally, we note that the gauge field A^ appears undifferentiated in 5v a = 
a -i jabcj^ x w h ereas it i s differentiated in v a = a~ 1 d ■ A a . This means that the 
perturbation SH-pp is softer than the unperturbed hamiltonian Hpp or, in other words, less 
divergent in the ultraviolet. If we calculate with the original hamiltonian, we get a certain 
number of divergent constants in the correlators. The result of a gauge transformation 
e a = rj^Xfj, on these correlators must agree with a calculation using the soft perturbation. 
This constrains the divergent constants. 
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Appendix D. Solution for and r( n ) 

The solution for and higher coefficient functions is similar to the solution for 
found in sec. 7. We differentiate (6.3) with respect to A x . four times and obtain, after 
setting A = 0, 

^xi,x (ri,x2,x3,x 4 ~l~ Qx;x2,x3,x 4 ) ^ ] part(4, 1) 

+ rg) xa>!B (ri 3 ) 3)X4 + Qg) 3 K j + ^ part (4, 2) = o, 

where 1?( 2 ) and are already known, and we have again used 1?( 2 ) = — Here 
part(n, n\) is the sum over all partitions of the set of n objects, x±, x 2 , ■ ■■x n , into subsets 
of n\ and n 2 = n — ri\ objects. In terms of the fourier transforms 

Q^CiZ!^^ 1 ^ 2 '^ 3 '^ 4 ) = (2tt)" 12 / d 4 k 1 d 4 k 3 d 4 k3d 4 k 4 exp(i^2k i ■ Xi) 

J i=i 

x 5(h + k 2 + h + k 4 ) Q^\z% a ; 4 (h-, k 2 , k 3 , k 4 ) 

(D.2) 

T(4) l\ a i? 2 a ^ 4 ( x i, x 2,X3,X4) = (27r)" 12 / d 4 k 1 d 4 k 3 d 4 k3d 4 k4exp(i^2k i ■ Xi) 



i=i 

x 5(h + k 2 + k 3 + k 4 ) fc 2, h, k 4 ), 

(D.3) 

where Q (3) £^ 3 3 £ 4 (fci; k ^ k ^ M and f {4) l\%f^S k ^ k ^ k ^ k ±) are defined onl y for k ^ + 
k 2 + ks + k 4 = 0, the equation for reads 

f £U*i) f (4) ^X 4 4 (fci, **, fcs, M + ^ P art(4, i) = - ff< 4 >-;££* (A*, fc, fca, 

(D.4) 

where 

H {4) Z%% a ^ **. h, = f g^fcx) Q (3) r 2 2 ; 3 3 ; 4 4 (^i; fc, *a, **) + J>art(4, l) 

+^ (4) rM 2 2 rA( fc i^2^3^ 4 ) 

x [ f (3 te(-*i " fc 2, h, h) + Q {2) Zl%(- k i ~ k; h, h) ] 

+ ^part(4,2). 

(D.6) 
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To solve (D.4), we project on each argument with a transverse transverse or longitu- 
dinal projector to obtain 



f (4)TTTT r; 2 r; 4 ( fc i^2,fe,fc 4 ) = - mi) +nkD +T(kD +T(ktr i 

f (4)LTTT ™; 4 4 ( fc i' k ^ fc 3, k 4 ) = - \a- l L(kl) + T(kl) + T(kl) + T(kj)]- 1 
etc. 

The formula for r( n ) for arbitrary n is similar. Each is expressed explicity and 
uniquely in terms of Q^ -1 -* and of ]?( 2 ) to r( n_1 ) which are already known. It is given 
by a symmetrized sum of products of two factors, as in eqs. (D.5) and (D.6), to which 
is applied a transverse or longitudinal projector onto each argument, and a division by 
E^rf^/cf), where vf\k 2 ) = T(k 2 ) or rf } (fc 2 ) = a" 1 L(A; 2 ). This gives all the 
uniquely in terms of to Q( n ~ r K 



Appendix E. Evaluation of T K t ' 

(3) 

We evaluate F K t , using the formulas of sec. 7, with the substitutions (9.14). 
From (7.15) and (8.6), we obtain 

iy S ) M2M 3 ( /c i' fc 2, fo) = a _1 f ^(/ci) ^t^^M^i; fc 2, fe) + (cyclic) 

= ia-'g ( (fc 3 ) M [ fgUfci) - (1 - 2) ] + (cyclic) ). ^ 

We apply transverse or longitudinal projectors to each Lorentz index, and use 
rg(fe) = T(P)P^(k) + a^L{k 2 )P^{k) to obtain 

^ (3) IJJ 3 ( fc 1^2,fe) = 0, 

^Ji^^) = ^ (T X -T 2 ) (fc 3 ) M 3 [P T (/C1)P T (/C 2 )] M1M2 , 

P^i^^s) = ia-'g ( (a-% -Ti) (/c 2 ) M2 [P T (fci)P L (fc 3 ) W 3 - (2 - 3)) 
H^l 1 l^ 5 {k 1 MM) = ia- 2 g ((L 2 -L 3 ) [P L (fc2)P L (fe 3 )] M w + cyclic), 

(E.2) 
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where we have used the notation Tj = T(k 2 ) and Li = L(k 2 ). From (7.16) and (7.17), we 
obtain finally 

rg ) IJJ3( /c i^2,/c 3 ) = o, 

rg^i.^, **,*&) = -ig aTi T +~ T ^ 2 +L3 [P T (h)p T (k 2 )]^ 

f g?^3( fc i^2,fc3) = - ^a- 1 ^ 7 ( Li L *^ L * L3 (*0mi [^(feO^faOWs + (cyclic)) . 



Appendix F. Evaluation of loop integrals 

We evaluate the integral that appears in (11.4) namely 



1 f, d , kjp-ih-kr 



J " (A;2)«t+2 ( 2yr )rf y - (k\y+^ [{k - ki) 2 ] l+olL 

We write this as 

Ir = m^w^m^T) r dadl3 a ° L flT ' (R2) 



where 



# T = (2iv)- d J d d k x [k 2 k 2 - (h ■ k) 2 ] exp[ - ak 2 - (3{k x - k) 2 ]. (F.3) 
We complete the square in the exponent, 

ak\ + f3(h - k) 2 = (a + f3)p 2 + (a + P^a^k 2 , 

where p = k\ — (a + (3)~ 1 f3k, and obtain 

R T = exp[ - (a + fi^afik 2 ] (27r)" d J d d p [p 2 k 2 - (p ■ k) 2 ] exp[ - (a + (3)p 2 ] 



(d-1) k 



2 



exp[ - (a + P) 1 a/3k 2 



2 (4n) d / 2 (a + /3) 1 + d / 2 

(FA) 

This gives 

lT = 2 (47r) d / 2 (P)^+i r 2 (l + a L ) ^ (F ' 5) 
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where 



S T = J dad(3 - ^ exp[ - (a + P^aPk 2 ]. (F.6) 

We insert the identity 1 = / °° da 5(a + f3 — cr), and change variables according to a = era' 
and P = a/3'. This gives, after dropping primes, 

/■OO 

S T = dad(3da 8(a + p- 1) a Qi p^ a 2 ^~d/ 2 exp [ _ a/3aA; 2 j 
Jo 

rOO 

= (/c 2 ) d/2_2ai_1 T(2a L + 1 - d/2) / dad/3 5(a + /3- 1) o^ 2 "^" 1 p*/ 2 -"*-- 1 

Jo 



= ^2^/2-20,-1 r(2a L + l-d/2) r 2 (rf/2- aL ) 
1 J r(rf-2a L ) 

We obtain finally, 

r(2a L + l-rf/2) T 2 (rf/2-a L ) 

i T = 



2 (47r) d / 2 T 2 (l + a L ) T(d - 2a L ) 
r(2a L -l) T 2 (2-« L ) 



(F.7) 



(F.8) 



"2(4tt) 2 T 2 (l + a L ) r(4-2a L )' 

for d = 4. 

We also evaluate the integral that appears in (11.9), 



- 2 (f.9) 



^ (k 2 ) 1+olL (2-K) d J 1 (k 2 ) 2 + a ^ [{h -k) 2 }^^ -k) 2 + k 2 ]' 

It contains the denominator [(k\ — k) 2 + k 2 } that comes from the non-local vertex. This 
integral is convergent in the ultraviolet for d < 4 + 2ar + 2a l. We shall evaluate it for d 
satisfying this condition, and then continue in d. We write it as 



- 2 



u (k 2 y+^r(2 + a T ) r(a L ] 

where 



p OO 

/ dadpd-f a 1+aT p aL ~ x eM-lk 2 ) R L , (F.10) 
Jo 



R L = (27i)- d J d d h [k\ k 2 - (ki ■ k) 2 ] exp[- akj - (/3 + ^(h - k) 2 ]. (F.ll) 
We complete the square in the exponent, 



ak\ + (P + 7 )(fci - k) 2 = (a + p + 7 )p 2 + (a + P + 7 ) _1 «(/3 + 7 )£; 2 , 
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This gives 



where p = k\ — (a + (3 + 7) 1 (/3 + 7 )fc, and obtain 

R L = exp[ - (a + /3 + 7) _1 a(/3 + 7) k 2 ] 

x (2n)- d J d d p [p*k 2 -(p- kf ] exp[ -(a + (3 + 7 )p 2 ] (pu2) 

2(4,)^+^ 7 )^/ 2 6XP[ " ( " + " + 7)_la(/3 + ] " 

L (P)^ (47r) d / 2 r(2 + a T )r(a L ) L ' 1 J 

/•oo l+a T flaL-l 

S L = / efc^ 7 -, ^n/, 1+d/2 exp[ - 7 fc 2 - (a + (3 + 7 )" 1 «(/3 + 7 )fc 2 ]• (F-14) 
Jo (« + P + 7J 1+C(// 

We insert the identity 1 = J °° da 5 (a + (3 + 7 — a) , and change variables according to 
a = era', (3 = a(3' and 7 = <7 7 '. This gives, after dropping primes, 

poo 

S L = dad(3d~fda S(a + (3 + 7 - 1) a 1+QT /3^-V +QT+Q ^- d/2 

x exp{ -a[ 7 + a(/3 + 7)]fc 2 } ( FJ5 ) 
= (k 2 )- aT - aL - 2+d/2 r(2 + a T + a L - d/2) M L , 

The argument of the T-function is positive in the region of convergence of the integral, 
d < 4 + 2ar + 2ctL- Here Ml is the finite integral 



poo 

= / dad/3d-f 8{a + f3 + 1 - 1) a 1+aT /3 ai_1 [a(/3 + 7) +7]°^ 
= d/3 da a 1+aT f3 aL ~ 1 {l -a 2 - f3) aL , 



Mi 



JO JO 

where we have used ax + 2a l = — (4 — d)/2. This gives 



{d-1) Fj-OL 

(4yr) d / 2 F(2 + a T )T(a L ) 



il = - :l , m l 



(d-l)T(l-a L ) Ml (FJ7) 



(47r) d / 2 r(2 + « T )r(l + a L ) 



Note that II is negative in the region of convergence of the integral, but after the contin- 
uation in <i, it is positive. 
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To evaluate M L we change variable to x = a 2 , and obtain 

<i-P) 2 



M, 



,1 Ai-py 
(1/2)/ d/? / ds;W 2 r^fl-/]-!)^, 

Jo Jo 



and upon changing variables to x = (1 — (3)y, we get 

.1 ,1-/3 



M L = (l/2) [ dpi dyy aT / 2 P aL - 1 (l-P) 1+aL+aT / 2 (l-y) aL . 
Jo Jo 

We again use «t + 2«l = — (4 — gQ/2 to write this as 

M L = (l/2) [ dy y aT/2 (l-y) aL [ * dp p^'^l - P) d / A . 
Jo Jo 

This is integrable by quadrature for d = 4, and in this case it gives 

»1 ,(l-y)2<*L (l-yf<*L+l 



M L = (1/2) J dy y~^ ( 
= (V2)( 



a L cxl + 1 

T(l-a L ) r(l + 2a L ) _ T(l-a L ) T(2 + 2a L ) 
a L T(2 + a L ) (a L + 1) T(3 + a L ) 

(-a 2 + 2a L + 2) T(l - a L ) T(2a L + 1) 



2a L (a L + l) T(a L + 3) 
where we used olt = — 2a^. This gives finally 

_ 3 (-a| + 2a L + 2) T 2 (1-« L ) T(2« L + 1) 



2 (4tt) 2 a L T(2 - 2a L )I> L + 2) F(a L + 3) ' 
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Figure Captions 



Fig. 1. Fig. 1. Diagrammatic representation of the functional DS equation for the quan- 
tum effective drift force, Q(A), in the presence of external sources A, eq. (6.1). 
The vertices are the tree-level vertices of the drift force K. The internal lines 
represent the exact gluon propagator T>(A) in the presence of the external source. 
The circle is the exact 3-gluon vertex of the quantum effective action T( A) in the 
presence of the external source. 

Fig. 2. Fig. 2. Diagrammatic representation of the DS equation for the gluon propagator, 
eq. (8.1). The vertices are the tree level vertices of the drift force K. The internal 
lines are the exact gluon propagator D with sources set to 0. The circles represent 
the exact 3 and 4-gluon vertices of the quantum effective action T, with sources 
set to 0. 
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Q( 1 )(x;y) = 
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+1/2 
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